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ABSTRACT 



A method for digital clock recovery and selective filtering 
includes prescribing or calculating first coefficients of a 
prototype of a selective filter at a characteristic frequency fc 
for a given sampling frequency fa. Second coefficients of a 
selective filter are calculated at the characteristic frequency 
fc2 from the first coefficients at points t k =At+k*d, where k 
is an integer greater than or equal to 0, by interpolating 
values of a continuous-time impulse response at points t k , 
where 



fc2*fa 



fc*fa { ' 



The selective filter is operated with the second coefficients 
at fa r 
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METHOD FOR SELECTIVE FILTERING 

CROSS-REFERENCE TO RELATED 
APPLICATION 

[0001] This application is a continuation of copending 
International Application No. PCT/DE99/01877, filed Jun. 
28, 1999, which designated the United States. 

BACKGROUND OF THE INVENTION 
10002] 1. Field of the Invention: 

l[0003] The invention relates to a method for selective 
filtering with a coefficient an d filte^ calculation unit, a phase 
detector, a loop filter, and a digital oscillator, these units 
simulating different selective filters depending on coeffi- 
cients used for the filter calculation. Systems for receiving 
digital data generally require a device for recovering the data 
clock signal implicitly contained in the signal. In principle, 
there are two approaches for receiving the data: (1) the clock 
control loop contains the circuit that supplies the sampling 
clock of the analog/digital converter ("ADC") (e.g. VCXO); 
and (2) the clock signal of the ADC is asynchronous with 
respect to the data clock. A fully digital contro l circuit 
calculates by interpolation from the output va±ae?bf the 
ADC samples in a timing pattern that is synchronous with 
the data clock. There exist prior art circuits that achieve the 
task in baseband. Most applications have filters upstream or 
downstream 61 the ADC, which filters select the signal that 
is intended to be-pfeg gssed, i.e., suppress adjacent signals 
and other interfeiejicfi_signals. These filters are generally 
optimized for a fixed bandwidth or data rate, such as, e.g., 
surface acoustic wave filters. Stringent requirements regard- 
ing phase linearity are usually imposed on the filter so that 
the pulse shape of the data signals is not distorted^ 

[0004] Basic theories with regard to fully digital clock 
recovery in digital modems by interpolation are discussed in 
an article by Floyd M. Gardener titled "Interpolation in 
Digital-Modems -Part I: Fundamentals". The article 
appeared in the journal IEEE Transactions on Communica- 
tions, Vol. 41, No. 3, in March 1993. 

[0005] A further approach in the prior art uses analog 
filters having a controllable bandwidth upstream of the 
ADC. However, these filters are expensive and, particularly 
in CMOS technology, are difficult to integrate on a circuit. 
Analog filters fundamentally have phase distortions that 
have to be reduced using additional circuits. Moreover, the 
sampling frequency of the ADC has to be adapted to the 
signal bandwidth. Another approach lies in providing a 
multistage digital selection filter upstream of the clock 
recovery and in carrying out the fine adjustment of the 
sampling frequency by interpolation again. After each stage 
of the filter, the sampling rate is reduced by a fixed factor. 
The approach has a disadvantage in that many filter stages 
are required for small bandwidths. Thus, the filter becomes 
complicated. Moreover, a fixed gradation means that not all 
the interference components can be suppressed, which gen- 
erally leads to an increase in outlay in the downstream 
interpolator. 

[0006] Lambrette U et al: "VARIABLE SAMPLE RATE 
DIGITAL FEEDBACK AND TIMING SYNCHRONIZA- 
TION" in I.E.E.E. Vehicular Technology Conference, New 
York, USA, I.E.E.E., Bd. Conf. 47, pages 1348-1352, dis- 



cusses two algorithms for digital receivers for processing a 
broader range of different sampling rates. One of the algo- 
rithms is also based on filtering the received signals prior to 
the time synchronization. A time synchronization algorithm 
is presented that is not data-aided, is based on digital 
feedback, and can process symbol rates deviating from a 
sampling rate. 

[0007] The paper by D. Kim et al.: "DESIGN OF OPTI- 
MAL INTERPOLATION FILTER FOR SYMBOL TIMING 
RECOVERY" in I.E.E.E. Transactions on Communications, 
I.E.E.E. Inc., New York, USA, Vol. 45, No. 7, pages 877- 
884, discloses an optimized interpolation filter for recover- 
ing the symbol timing in a digital receiver, in which the 
sampling rate of the analog-to-digital converter on the input 
side is not synchronized with the symbol clock of the 
transmitter. 

[0008] The paper by K Bucket et al. "PERIODIC TIM- 
ING ERROR COMPONENTS IN FEEDBACK SYN- 
CHRONIZERS OPERATING ON NONSYNCHRONIZED 
SIGNAL SAMPLES" I.E.E.E. Transactions on Communi- 
cations, I.E.E.E. Inc., New York, USA, Vol. 46, No. 6, pages 
747-749, reveals that the synchronization error contains 
periodic components through a loop for timing recovery on 
detection of nonsynchronized samples of a noisy sine signal. 
These periodic errors are produced exclusively by non-ideal 
interpolation between the nonsynchronized signal samples 
and disappear when synchronized sampling is performed. 

SUMMARY OF THE INVENTION 

[0009] It is accordingly an object of the invention to 
provide a method for selective filtering that overcomes the 
hereinafore-mentioned disadvantages of the heretofore- 
known devices and methods of this general type and that can 
be realized with little outlay. 

[0010] With the foregoing and other objects in view, there 
is provided, in accordance with the invention, a method for 
selective filtering, including the steps of simulating different 
selective filters with a coefficient and filter calculation unit, 
a phase detector, a loop filter, and a digital oscillator 
depending on first and second coefficients used for filter 
calculation, determining the first coefficients of a prototype 
of a further selective filter at a characteristic frequency fc 
for — given sampling frequency fa, calculating the second 
coefficients of the further selective filter at a characteristic 
frequency fc2 from the first coefficients at points t k =At+k*d, 
where k-0, 1, . . . , the further selective filter to be operated 
at a sampling rate fa a , by interpolating values of a continu- 
ous-time impulse response of a simulation of the further 
selective filter at points t k , where 



= fc2*Ja 
" fc*f<*i ' 



[0011] simulating the further selective filter with the sec- 
ond coefficients using the coefficient and filter calculation 
unit, the phase detector, the loop filter, and the digital 
oscillator, and operating the further selective filter at the 
sampling rate fa 3 . 

[0012] The approach employed by the invention is to use 
a filter disclosed in International PCT publication WO-A- 
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00/02311, filed Jan. 13, 2000, corresponding to PCT/DE 
99/01878 and U.S. patent application having U.S. attorney 
Docket No. S0439-SB, by the same inventor and to combine 
it with a variant of a prior art clock recovery circuit with an 
interpolation filter. As such, the approach of deriving filter 
coefficients from the coefficients of a prototype filter is 
extended in order to additionally realize a variable delay that 
is required for the interpolation of the data signal and that 
must generally be shorter than the period of the ADC clock 
signal. Moreover, use is made of the fact that the filter 
bandwidth is always proportional to the data rate, i.e., the 
bandwidth and the decimation factor of the filter are in a 
fixed relationship with respect to one another. Correspond- 
ingly, in qualitative terms, there is more time available for 
calculating an output value of the filter, the lower the data 
rate or the bandwidth of the filter. On the other hand, it is 
also the case that more coefficients are required for calcu- 
lating the output value given a smaller bandwidth. Meaning, 
in quantitative terms, that the same number of arithmetic 
operations per second is always required regardless of the 
data rate or the filter bandwidth. A suitable circuit is 
described below. A variant that is optimized in respect of 
outlay is likewise presented, for a case where the operating 
clock of the circuit is higher than the sampling clock at the 
- filter input. 

[0013] The circuits described below accomplish two tasks: 

(1) to realize low-pass filters having a bandwidth propor- 
tional to the data rate for suppressing adjacent signals; and 

(2) generating an output signal that is phase-synchronous 
with the data clock by interpolation because the sampling 
clock of the input signal is not synchronous with the data 
clock (free-running oscillator). 

[0014] In order to accomplish the second task, a control 
loop including a phase detector, a loop filter, and a digital 
oscillator ("DTO") is used. The oscillator is realized as an 
overflowing accumulator. If an overflow occurs, then a 
sample is interpolated from the samples of the input signal 
and feeds the phase detector and downstream circuits. In the 
steady-state condition, some of these interpolated values 
lie — as desired — in the center of the data pulses. The inter- 
polated signal usually has a sampling frequency correspond- 
ing to twice the data rate (or symbol rate in the case of 
two-dimensional data transmission such as, e.g., QPSK or 
QAM). Because the two clock signals are asynchronous, the 
sampling instants of the signal to be interpolated generally 
lie between the sampling instants of the input signal. There- 
fore, a signal is derived from the state of the DTO after the 
overflow, the derived signal specifying the time interval 
between the desired sampling instant and the last sampling 
instant of the input signal. The input signal must be tempo- 
rally shifted by the interval amount by the interpolation 
filter. 

[0015] The invention accomplishes the second task by 
taking into account, during the interpolation of the filter 
coefficients from the coefficients of a prototype, not only the 
desired bandwidth but also the desired temporal shift. This 
. also indicates how the invention accomplishes the first task. 
The bandwidth of the decimating interpolation filter must be 
adapted to the data rate. For example, the output signal of the 
loop filter can be used for such a purpose. In the steady-state 
condition, the output signal is proportional to the difference 
between the desired interpolation rate and the idle frequency 
of the DTO. Because the idle frequency is prescribed, it is 



possible to form a signal that is exactly proportional to the 
desired bandwidth of the filter. The loop filter has an integral 
element whose output signal has, in principle, the same 
properties, coupled with the additional advantage that the 
signal has significantly less noise than the output signal of 
the overall filter. 

[0016] If, for a bandwidth fcl of the decimating interpo- 
lation filter, Ml coefficients are required for calculating 
output values at the rate fal, then M2-Mlfcl/fc2 coeffi- 
cients are required in the proposed method for calculating 
the filter coefficients for a bandwidth fc2 and an output rate 
fa2=fal*fc2/fcl. It follows, therefore, that the same pro- 
cessing speed is required in both cases, because fal*Mlo 
fa2*M2. Correspondingly, it is possible, in principle, always 
to carry out the filter calculation with the same number of 
arithmetic elements, regardless of the data rate. The proto- 
type is dimensioned for a bandwidth fc at a sampling 
frequency fa. In principle, the two characteristic frequencies 
are arbitrary, but in individual cases fc must be chosen to be 
low enough that, taking account of the chosen method for 
interpolating the coefficients, the accuracy requirements of 
the respective application are met. If the decimating inter- 
polation filter is intended to have a bandwidth fc2 — match- 
ing the sampling rate fa2 after decimation — and to be 
operated at a sampling frequency fal, then an "expansion 
factor" d=(fc2/fc)*(fa/fal) results, i.e., the sampling pattern 
of the prototype must be "expanded" by the factor d. 

[0017] If the intention is to interpolate the samples of the 
data signal at the sampling frequency fa2, then the DTO of 
the control loop must be operated with the increment 
AI=fa2/fal (AI is supplied by the DTO, see above). As such, 
it is assumed that the DTO is operated at the sampling 
frequency fal and overflows at the value 1 (i.e., the output 
values of the DTO He between 0 and 1). The output value 10 
after the overflow can have the value AI at most. It is defined 
that the coefficients of the prototypes are situated at points 
t=n where n=0, 1, etc. Correspondingly, the coefficients of 
the required filter that are to be interpolated are situated at 
the points t-At+k*d, where k-0, 1, etc., and the value At is 
defined by the control loop. The term k*d can be converted 
into k*AI*(fa/fc)*(fc2/fa2)=k*AI*r, where r is a constant 
factor (fa/fc is defined in the configuration of the prototype 
and problems dictate that fc2/fa2 is a constant). It likewise 
follows that t»At+k*d=(I0+k*AI)*r, where I0+k*AI— disre- 
garding noise terms — represent the output values of the 
DTO in the steady-state condition. It is thus the case that 
At«IO*r, in particular, holds true, i.e., the output value of the 
DTO after an overflow is multiplied by the constant r and is 
then used as a start value for the interpolation of the filter 
coefficients. 

[0018] In accordance with another mode of the invention, 
the determining step is performed by prescribing or calcu- 
lating first coefficients of a prototype of a further selective 
filter at a characteristic frequency fc for a given sampling 
frequency fa. 

[0019] In accordance with a further mode of the invention, 
output values I0+k*AI are produced with the digital oscil- 
lator, and using the output values I0+k*AI for calculating 
At+k*d. 

[0020] In accordance with an added mode of the inven- 
tion, a sampled data signal is selectively filtered and an 
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output value of the digital oscillator after an overflow 10, 
multiplied by the constant 



famfc2 
r= fc*fa2* 



[0021] is used for a start value At for an interpolation of 
the second coefficients (h2(t)), where fa2 represents a sam- 
pling frequency at which the data signal is interpolated. 

[0022] In accordance with an additional mode of the 
invention, an overflowing accumulator is used as the digital 
oscillator. 

[0023] In accordance with yet another mode of the inven- 
tion, the accumulator is fed with an increment formed by 
adding an output signal of the loop filter to a value (IncrO) 
defining an idle frequency of the digital oscillator. 

[0024] In accordance with yet a further mode of the 
invention, an input signal is fed to arithmetic units operating 
in parallel, each of the arithmetic units calculating a future 
output value of a selective filter, and an output value of an 
arithmetic unit having already finished a calculation at the 
time of the overflow is selected in the event of an overflow 
of the digital oscillator. 

[0025] In accordance with yet an added mode of the 
invention, a number N of arithmetic units operating in 
parallel is defined based on a number K of coefficients of the 
selective filter simulated respectively by the coefficient and 
filter calculation unit, the phase detector, the loop filter, and 
the digital oscillator, according to the following conditions: 

fafc2 



[0026] where N is an integer. 

[0027] In accordance with yet an additional mode of the 
invention, the arithmetic units are selected with a modulo-N 
counter, the arithmetic units are cyclically addressed with 
the modulo-N counter, and the modulo-N counter is incre- 
mented with each overflow of the digital oscillator. 

[0028] In accordance with again another mode of the 
invention, a sampled data signal is selectively filtered, an 
accumulator in each arithmetic unit is reset if the arithmetic 
unit is addressed, the arithmetic unit is loaded with an output 
value of the digital oscillator after the overflow, an instan- 
taneous increment value of the digital oscillator is stored, the 
instantaneous increment value is accumulated for each sam- 
pling clock of the input signal, an output value of the 
accumulator is scaled with a value T,131 for each input 
clock, where fa2 represents a sampling frequency at which 
the data signal is interpolated. 

[0029] In accordance with again a further mode of the 
invention, a position is defined with a scaled output value of 
the accumulator at which a second coefficient must be 
interpolated for the further selective filter from the first 
coefficients of the prototype. 

[0030] In accordance with again an added mode of the 
invention, a single coefficient calculation unit is provided, 



the second coefficients are calculated sequentially in time 
with the single coefficient calculation unit, the filter output 
values are calculated sequentially in time in a multiplier and 
splitting the filter output values between N accumulators, 
and the respective input values are selectively stored in N 
further accumulators. 

[0031] In accordance with a concomitant mode of the 
invention, the respective input values are an increment and 
output value of the digital oscillator. 

[0032] Other features that are considered as characteristic 
for the invention are set forth in the appended claims. 

[0033] Although the invention is illustrated and described 
herein as embodied in a method for selective filtering, it is, 
nevertheless, not intended to be limited to the details shown 
because various modifications and structural changes may 
be made therein without departing from the spirit of the 
invention and within the scope and range of equivalents of 
the claims. 

[0034] The construction and method of operation of the 
invention, however, together with additional objects and 
advantages thereof, will be best understood from the fol- 
lowing description of specific embodiments when read in 
connection with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWING 

[0035] FIG. 1 is a block circuit diagram for carrier recov- 
ery with a decimating interpolation filter according to the 
invention; 

[0036] FIG. 2 is a block and schematic circuit diagram of 
a digital values of FIG. 1; 

[0037] FIG. 3 is a block circuit diagram of a combined 
unit for calculating both the filter coefficients and the filter 
output values of FIG. 1; 

[0038] FIG. 4 is a block circuit diagram of an arithmetic 
unit from FIG. 3; 

[0039] FIG. 5 is a block circuit diagram of a MAC unit 
from FIG. 4; and 

[0040] FIG. 6 is a block circuit diagram of a modified 
arithmetic unit for the case where at least N operating clocks 
are available per clock of the input signal. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

[0041] In all the figures of the drawing, sub-features and 
integral parts that correspond to one another bear the same 
reference symbol in each case. 

[0042] Referring now to the figures of the drawings in 
detail and first, particularly to FIG. 1 thereof, there is shown 
a basic circuit diagram for carrier recovery with a decimat- 
ing interpolation filter. The control loop for clock recovery 
is illustrated including a coefficient and filter calculation unit 
10, a phase detector 12, a loop filter 14, and a digital 
oscillator 16. In the configuration, the input signals are fed 
to the coefficient and filter calculation unit 10, which cal- 
culates the output signals from the input signals. The output 
signals are fed not only to the output but also to the phase 
detector 12, which is, in turn, connected to the loop filter 14. 
The output of the loop filter 14 is fed to the digital oscillator 
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16. The digital oscillator 16 outputs to the coefficient and 
filter calculation unit 10 the signals "DTO output", "over- 
flow indication," and "increment" for controlling the calcu- 
lation unit 10. 

[0043] FIG. 2 further illustrates the digital oscillator 16 
(DTO). At its heart, the DTO 16 includes an accumulator 20 
with overflow ("wrap around", sawtooth-waveform output 
signal). The accumulator 20 is fed by an increment formed 
by adding the output signal of the loop filter 14 to a value 
IncrO. The value IncrO defines the idle frequency of the DTO 
16. The DTO increment is passed on to the coefficient 
calculation unit 10. (The value IncrO+integral element of the 
output signal of the loop filter can also be output here.) 
Furthermore, in the case of each accumulator overflow, a 
new calculation of a filter output value is initiated by the 
overflow indication 22 of the accumulator 20. 

[0044] Moreover, the DTO 16 includes a summing unit 24, 
in which the output of the loop filter 14 and the value IncrO 
are added. The value "increment", which corresponds to the 
current sum of the output of the loop filter 14 and the value 
IncrO, can be picked off at the output of the adder 24. The 
content of the accumulator 20 is available at the DTO output 
26. 

[0045] FIG. 3 details the coefficient and filter calculation 
unit 10, which serves for calculating the filter output values. 
The input signal feeds a plurality of arithmetic units 30 
operating in parallel. Each of these units 30 calculates a 
future output value of the filter. The number N of units 30 
required is defined by the number of coefficients of the 
prototype filter (N^l/r* (number of coefficients of the 
prototype filter) and N is an integer. If a DTO overflow is 
signaled, the output value of the arithmetic unit 30 that has 
already finished the filter calculation at this point in time is 
copied into the result memory 34 through a selector switch 
32. The arithmetic unit 30 is subsequently initialized. At its 
heart, the control unit 36 has a "modulo counter," which 
cyclically addresses the arithmetic units 30. The modulo 
counter is incremented with each DTO overflow. 

[0046] Feeding the input signal to all N arithmetic units 30 
configures the coefficient and filter calculation unit 10. The 
output value 26 of the DTO 16 and the increment signal 
thereof are likewise fed to all of the units 30. Also provided 
is a control unit 36 with a modulo counter and N outputs, 
each of which is selectively connected to one of the arith- 
metic units 30. The modulo counter in the control unit 36 is 
advanced with the overflow indication of the DTO 16. At the 
same time, in the case of each overflow, the command "load" 
is output for the result memory. 

[0047] FIG. 4 illustrates the arithmetic unit 30 in detail. 
The arithmetic unit 30 includes a MAC unit 40 that is 
illustrated in more detail in FIG. 5. The MAC unit 40 is fed 
with the input signal and, from this, supplies an output signal 
to the selector switch 32. Furthermore, the arithmetic unit 30 
has an accumulator with a scaler that buffer-stores and scales 
the values Increment of the DTO 16 and output value 26 of 
the DTO 16. These values are then fed to a coefficient 
calculation unit 44, which supplies the coefficients for the 
MAC unit 40. For such purpose, the coefficient calculation 
unit 44 is connected to a coefficient memory 46, from which 
the corresponding coefficient values for the prototype filter 
can be read out respectively. 

[0048] The construction of the MAC unit 40 is illustrated 
in detail in FIG. 5. The MAC unit 40 includes a multiplier 



52, to which the input signal and the corresponding coeffi- 
cient from the coefficient calculation unit 44 are fed. The 
output value of the multiplier is stored in the accumulator 
50. As a result, the accumulator 50 makes the output signal 
available permanently, which signal is then fed to the 
selector switch 32. The accumulator 50 is reset by the 
selection/initialization command from the control unit 36. 

[0049] The arithmetic unit 30 functions as follows. If the 
arithmetic unit 30 is selected by the control unit 36, (1) the 
accumulator 50 of the MAC unit 40 is reset, (2) the 
accumulator 42 upstream of the coefficient calculation unit 
44 is loaded with the instantaneous output value 26 of the 
DTO 16 (value after the overflow), and (3) the instantaneous 
increment value of the DTO 16 is stored. After the initial- 
ization, the increment is accumulated in the case of each 
sampling clock of the input signal. 

[0050] The output value of the accumulator 42 is scaled 
with the value "r" in the case of each input clock. The result 
defines the position at which a coefficient must be interpo- 
lated from the coefficients of the prototype filter (the inter- 
polation method is arbitrary, in principle). The calculated 
coefficient is multiplied by the current filter input value and 
is subsequently accumulated. 

[0051] FIG. 6 shows a modified embodiment of the inven- 
tion that can be used when at least N operating clocks are 
available per clock of the input signal. In such a case, N 
arithmetic units 30 are not required. Rather a single multi- 
plier 152 and a single unit 144 for calculating the filter 
coefficients (including ROM 146 and scaler) suffices. 

[0052] Thus, FIG. 6 shows a modified coefficient and 
filter calculation unit 10, which, just like the coefficient and 
filter calculation unit 10 described above in FIGS. 2 to 5, can 
be used at the corresponding point in FIG. 1. For the 
configuration, the input signal is fed only to a single mul- 
tiplier 152, but, from there, it is selectively split between N 
accumulators 150 and fed from the accumulators 150, in 
turn, to the result memory 134. The result memory 134 then 
always contains the output signal. A control unit 136 is also 
provided and is connected to the overflow indication of the 
DTO 16. The control unit 136 includes a "modulo counter" 
for driving the N accumulators 150. In this case, too, the 
multiplier 152 is driven by the coefficient calculation unit 
144, which once again reads the values of the prototype 
coefficients from the ROM memory 146. N accumulators 
142 are additionally provided for storing N different values 
for the increment of the DTO 16 and the output value 26 of 
the DTO 16. These accumulators 142 also receive their load 
command selectively from the control unit 136. 

[0053] One of the N accumulators drives the unit cycli- 
cally. The N accumulators downstream of the multiplier can 
be combined to form a circuit containing an adder and N 
memory cells. 



I claim: 

1. A method for selective filtering, which comprises: 

simulating different selective filters with a coefficient and 
filter calculation unit, a phase detector, a loop filter, and 
a digital oscillator depending on first and second coef- 
ficients used for filter calculation; 
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determining first coefficients of a prototype of a further 
selective filter at a characteristic frequency fc for a 
given sampling frequency fa; 

calculating second coefficients of the further selective 
filter at a characteristic frequency fc2 from the first 
coefficients at points t k =At+k*d, where k=0, 1, . . . the 
further selective filter to be operated at a sampling rate 
fa a , by interpolating values of a continuous-time 
impulse response of a simulation of the further selec- 
tive filter at points t k , where 



_ fc2*fa 
" fc*fai ' 



simulating the further selective filter with the second 
coefficients using the coefficient and filter calculation 
unit, the phase detector, the loop filter, and the digital 
oscillator; and operating the further selective filter at 
the sampling rate fa 2 . 

2. The method according to claim 1, which further com- 
prises performing the determining step by prescribing first 
coefficients of a prototype of a further selective filter at a 
characteristic frequency fc for a given sampling frequency 
fa. 

3. The method according to claim 1, which further com- 
prises performing the determining step by calculating first 
coefficients of a prototype of a further selective filter at a 
characteristic frequency fc for a given sampling frequency 
fa. 

4. The method according to claim 1, which further com- 
prises producing output values I0+k*AI with the digital 
oscillator, and using the output values I0+K*AI for calcu- 
lating At+k*d. 

5. The method according to claim 1, which further com- 
prises: selectively filtering a sampled data signal; and using 
an output value of the digital oscillator after an overflow 10, 
multiplied by the constant r=fa*fc2 /fc*fa2, for a start value 
At for a n interpolation of the second coefficients (h2(t)), 
where fa2 represents a sampling frequency at which the data 
signal is interpolated. 

6. The method according to claim 1, which further com- 
prises using an overflowing accumulator as the digital 
oscillator. 

7. The method according to claim 6, which further com- 
prises feeding the accumulator with an increment formed by 
adding an output signal of the loop filter to a value (IncrO) 
defining an idle frequency of the digital oscillator. 

8. The method according to claim 1, which further com- 
prises: 

feeding an input signal to arithmetic units operating in 
parallel, each of the arithmetic units calculating a future 
output value of a selective filter; and 

selecting, in the event of an overflow of the digital 
oscillator, an output value of an arithmetic unit having 
already finished a calculation at the time of the over- 
flow. 

9. The method according to claim 8, which further com- 
prises defining a number N of arithmetic units operating in 
parallel based on a number K of coefficients of the selective 
filter simulated respectively by the coefficient and filter 



calculation unit, the phase detector, the loop filter, and the 
digital oscillator, according to the following conditions: 



K-fc-fa2 
fafc2 

where N is an integer. 

10. The method according to claim 8, which further 
comprises: 

selecting the arithmetic units with a modulo-N counter; 

cyclically addressing the arithmetic units with the 
modulo-N counter, and 

incrementing the modulo-N counter with each overflow of 
the digital oscillator. 

11. The method according to claim 9, which further 
comprises: 

selecting the arithmetic units with a modulo-N counter; 

cyclically addressing the arithmetic units with the 
modulo-N counter; and 

incrementing the modulo-N counter with each overflow of 
the digital oscillator. 

12. The method according to claim 10, which further 
comprises: 

selectively filtering a sampled data signal; 

resetting an accumulator in each arithmetic unit if the 
arithmetic unit is addressed; 

loading the arithmetic unit with an output value of the 
digital oscillator after the overflow; 

storing an instantaneous increment value of the digital 
oscillator; 

accumulating the instantaneous increment value for each 
sampling clock of the input signal; 

scaling an output value of the accumulator with a value 



fc*fa2 



for each input clock, where fa2 represents a sampling 
frequency at which the data signal is interpolated. 

13, The method according to claim 11, which further 
comprises: 

selectively filtering a sampled data signal; 

resetting an accumulator in each arithmetic unit if the 
arithmetic unit is addressed; 

loading the arithmetic unit with an output value of the 
digital oscillator after the overflow; 

storing an instantaneous increment value of the digital 
oscillator; 
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accumulating the instantaneous increment value for each 
sampling clock of the input signal; 

scaling an output value of the accumulator with a value 



fa*fc2 
r= fe*fd2 



for each input clock, where fa2 represents a sampling 
frequency at which the data signal is interpolated. 

14. The method according to claim 12, which further 
comprises defining, with a scaled output value of the accu- 
mulator, a position at which a second coefficient must be 
interpolated for the further selective filter from the first 
coefficients of the prototype. 

15. The method according to claim 13, which further 
comprises defining, with a scaled output value of the accu- 



mulator, a position at which a second coefficient must be 
interpolated for the further selective filter from the first 
coefficients of the prototype. 

16. The method according to claim 1, which further 
comprises: 

providing a single coefficient calculation unit; 

calculating the second coefficients sequentially in time 
with the single coefficient calculation unit; 

calculating the filter output values sequentially in time in 
a multiplier and splitting the filter output values 
between N accumulators; and 

selectively storing the respective input values in N further 
accumulators. 

17. The method according to claim 16, wherein the 
respective input values are an increment and output value of 
the digital oscillator. 

* * * * * 



3/24/05, EAST Version: 2.0.1.4 



A 



^(SrtJnited States Patent 

( Petrus ) 



US006177906B1 

(io) Patent No.: US 6,177,906 Bl 
(45) Date of Patent: Jan. 23, 2001 



(54) MULTIMODE ITERATIVE ADAPTIVE 

SMART ANTENNA PROCESSING METHOD 
AND APPARATUS 

(75) Inventor: Paul Petrus, Sunnyvale, CA (US) 

(73) Assignee: ArrayComm, Inc., San Jose, CA (US) 

( * ) Notice: Under 35 U.S.C 154(b), the term of this 
patent shall be extended for 0 days. 

(21) Appl. No.: 09/286,135 

(22) Filed: Apr. 1, 1999 

(51) Int. CI. 7 G01S 3/16 

(52) U.S. CI 342/378; 342/373; 370/310; 

370/347 

(58) Field of Search 342/378, 380, 

342/382, 383, 367, 373; 370/310, 321, 
334, 347 

(56) References Cited 

U.S. PATENT DOCUMENTS 



5,515,378 * 5/1996 Roy et al 370/95.1 

5,887,038 3/1999 Golden 375/347 

5,909,470 * 6/1999 Barratt et al 375/324 

5,930,243 ♦ 7/1999 Parish et al 370/334 

6,018,317 * 1/2000 Dogan et al 342/378 



OTHER PUBLICATIONS 

Rias Muhamed and T.S. Rappaport: "Direction of arrival 
estimation using antenna arrays," Technical Report 
MPRG-TR-96-03, Mobile and Portable Radio Research 
Group, Bradley Department of Electrical engineering, Vir- 
ginia Polytechnic Institute, Jan. 1996, Section 3.8, pp. 64-78 
and Bibliography, pp. 151-158 (also available as Rias 
Muhamed: "Direction of arrival estimation using antenna 
arrays," Masters Thesis, Bradley Department of Electrical 
Engineering, Virginia Polytechnic Institute and State Uni- 
versity, Blacksburg, VA 24061, USA). 



B. G. Agee: "Blind separation and capture of communication 
signals using a multitarget constant modulus beamformer," 
1989 IEEE Military Communications Conference ("Milcom 
89"), vol. 2, pp. 340-346, New York: IEEE, 1989. 

J. Lundell and B. Widrow: "application of constant modulus 
adaptive beamformer to constant and non-constant modulus 
signals," Proceedings, 1988 Asilomar Conference on Sig- 
nals, Systems and Computers (ACSSC-1988), pp. 432-436, 
1988. 

C. B. Papadias and A. Paulraj, a space-time constant modu- 
lus algorithm for SDMA systems, Proceedings, IEEE 46th 
Vehicular Technology Conference, pp. 86-90, 1996. 

* cited by examiner 

Primary Examiner — Thomas H. Tarcza 

Assistant Examiner — Dao L. Phan 

(74) Attorney, Agent, or Firm — Dov Rosenfeld; Inventek 

(57) ABSTRACT 

A method and apparatus is described for implementing 
adaptive smart antenna processing in a receiving communi- 
cation station that includes an array of antennas and means 
for adaptive smart antenna processing, the method and 
apparatus including determining weight vectors for the 
adaptive smart antenna processing. Using the invention 
offers advantages when operating in a low SINR 
environment, for example, in a mobile environment in which 
the remote users are travelling at high speeds, hence the 
signals undergo fading. One aspect is hybrid weight adap- 
tation that starts off with a method with good convergence 
properties, for example, one known to converge in a low 
SINR environment then switches to a method that converges 
rapidly, for example when started with relatively high qual- 
ity initial conditions. To deal with high mobility, the weights 
determined from data at a particular burst are applied on that 
particular burst. Such weights may not be optimal for the 
subsequent bursts. When several users are present in a given 
channel, a multiport architecture is used to track each 
individual remote user. 
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MULTIMODE ITERATIVE ADAPTIVE 
SMART ANTENNA PROCESSING METHOD 
AND APPARATUS 

FIELD OF THE INVENTION 

This invention is related to wireless communication 
systems, and in particular to determining weights for adap- 
tive smart antenna processing in a wireless communication 
receiver having an array of antenna elements and means for 
adaptive smart antenna processing. 

BACKGROUND 

Wireless communications systems that include commu- 
nication stations that have an antenna array and means for 
adaptive smart antenna processing are known. Such com- 
munications stations are sometimes called smart antenna 
communications stations. When receiving a signal from a 
subscriber unit, the signals received by each of the antenna 
array elements are combined by the adaptive smart antenna 
processing means to provide an estimate of a signal received 
from a particular subscriber unit. With the smart antenna 
processing comprising linear spatial processing, each of the 
complex-valued (i.e., including in-phase I and quadrature Q 
components) signals received from the antenna elements is 
weighted in amplitude a nd phase by a weighting factor and 
the weighted sigaals are then summed to provide the esti- 
mate. The adaptive smart antenna processing means can then 
be described by a set of complex valued weights, one for 
each antenna elements. These complex valued weights in 
turn can be described as a single complex valued vector of 
m elements, where m is the number of antenna elements. 
This can be extended to include spatio-temporal processing, 
where the signal at each antenna element, rather than being 
simply weighted i n ampliti irfc an ^ p h r^ is filtered by some 
complex valued filter, typically for time equalization. Each 
filter can be described by a complex-valued transfer function 
or convolving function. The adaptive smart antenna pro- 
cessing of all elements can then be described by a complex 
valued m-vector of m complex valued convolving functions. 

Several methods are known for determining the weight 
vectors of received signals. These include methods that 
determine the directions of arrival of signals from subscriber 
units, and methods that use the spatial characteristics of 
subscriber units, for example, the spatial signatures. See for 
example U.S. Pat. Nos. 5,515,378 and 5,642,353 entitled 
SPATIAL DIVISION MULTIPLE ACCESS WIRELESS 
COMMUNICATION SYSTEMS, to Roy et al., for methods 
that use directions of arrival, and U.S. Pat. No. 5,592,490 
entitled SPECTRALLY EFFICIENT HIGH CAPACITY 
WIRELESS COMMUNICATION SYSTEMS, to Barratt et 
al„ and U.S. Pat. No. 5,828,658 entitled SPECTRALLY 
EFFICIENT HIGH CAPACITY WIRELESS COMMUNI- 
CATION SYSTEMS WITH SPATIO-TEMPORAL 
PROCESSING, to Ottersten et al., for methods that use 
spatial signatures. So-called "blind" methods determine the 
weights from the signals themselves, but without resorting 
to training signals — that is without determining what 
weights can best estimate a known symbol sequence. Such 
methods usually use some known characteristic of the signal 
transmitted by the subscriber unit to determine the best 
weights to use by constraining the estimate to have this 
property, and hence are called property restoral methods. 
Property restoral methods in turn can be classified into two 
groups. "Partial" property restoral methods restore one or 
more typically simple properties of the signal without com- 
pletely reconstructing the modulated received signal, for 
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example by demodulating and then re- modulating. "Deci- 
sion directed" (DD) methods construct an accurate copy of 
the signal by making symbol decisions (e.g., demodulating) 
the received signal. 

5 One example of the first group, partial restoral methods, 
is the constant modulus (CM) method, which is applicable 
to communications systems that use a modulation scheme 
that has a constant modulus, including, for example phase 
modulation (PM), frequency modulation (FM), phase shift 

10 keying (PSK) and frequency shift keying (FSK). See for 
example J. R. Treichler; M. L. Larimore: "New Processing 
Techniques Based on the Constant Modulus Algorithm," 
IEEE Transactions on Acoustics, Speech, and Signal 
Processing, vol. ASSP-33, No. 2, pp. 420-431, April 1985. 

j 5 Other partial property restoral techniques include techniques 
that restore the spectral properties of the signal, such as the 
spectral self-coherence. Spectral coherence restoral tech- 
niques use known spectral coherence properties of any 
signals received at the antenna array. For example, in certain 

20 situations, the signals may be assumed to be cyclo- 
stationarity, i.e., to have periodic autocorrelation functions. 
Other methods include those that restore high order 
statistics, e.g., moments or cumulants. See for example B. 
Agee, S. Schell, W. Gardner: "Spectral Self-Coherence 

25 Restoral: A New Approach to Blind Adaptive Signal Extrac- 
tion Using Antenna Arrays," Proceedings of the IEEE, vol. 
78, No. 4, April 1990, and U.S. Pat. No. 5,260,968 to 
Gardner, et al., entitled METHOD AND APPARATUS FOR 
MULTIPLEXING COMMUNICATIONS SIGNALS 

30 THROUGH BLIND ADAPTIVE SPATIAL FILTERING, 
and U.S. Pat. No. 5,255,210 to Gardner, et al., entitled 
SELF-COHERENCE RESTORING SIGNAL EXTRAC- 
TION APPARATUS AND METHOD. 

Decision directed methods use the fact that the modula- 

35 tion scheme of the transmitted subscriber unit signal is 
known, and determine weights that produce a signal (a 
"reference signal") that has the required modulation scheme, 
and if transmitted by a remote user, would produce signals 
at the antenna elements of the array that are "close" to the 

40 signals actually received, the reference signal production 
including making symbol decisions. See for example U.S. 
patent applications Ser. No. 08/729,390 entitled METHOD 
& APPARATUS FOR DECISION DIRECTED DEMODU- 
LATION USING ANTENNA ARRAYS & SPATIAL PRO- 

45 CESSING to Barratt, et al. (filed Oct. 11, 1996), and Ser. No. 
09/153,110 entitled METHOD FOR REFERENCE SIG- 
NAL GENERATION IN THE PRESENCE OF FRE- 
QUENCY OFFSETS IN A COMMUNICATIONS STA- 
TION WITH SPATIAL PROCESSING to Petrus, et al. (filed 

50 Sep. 15, 1998), for descriptions of systems that use decision 
directed weight determination. 

Some iterative methods, including partial restoral 
methods, for example, the CM method, are known to con- 
verge even for low signal-to-noise ratios (SNRs), low signal- 

55 to-interference-plus-noise-ratios (SINRs), and high fading 
situations as would be encountered in communication sys- 
tems wherein the subscriber units are highly mobile. Such 
methods are called "iterative weight determining methods 
with good convergence properties" herein. Methods with 

60 good convergence properties may however take many itera- 
tions to converge. The CM method, for example, may so 
take many iterations to converge, and therefore may not 
converge fast enough in an actual system. For example, in a 
high mobility system, it is desired to use the weight vector 

65 on a current burst that is derived from the current burst's 
data. This implies rapid calculation of the weights, which 
may not be possible with the CM method. The decision 
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directed method, on the other hand, is one example of a class have good convergence properties with the advantages of 

of methods that converges rapidly if the initial condition, for methods that converge rapidly. 

example, the initial signal-to-noise ratio (SNR), and signal- Another object of the present invention is a method and 
to-interference-plus-noise-ratio (SINR) is high, or the initial apparatus for "blind" weight determination that perform 
weight vector is sufficiently close to the correct value. 5 we u under low SINR and high fading situations, and that 
Methods that so converge rapidly if the initial weight vector converge rapidly, i.e., in a small number of iterations, 
is sufficiently close to the correct value are called "rapidly Qbject ^ fl ad ^ ^ ^ 
converging iterative weight determining methods herein. method and tus that determines adaptive smart 
Rapidly converging methoos such as the DD method are ^ wei hts efficieD ii y in a low signal-to- 
becoming more widely used in smart-antenna based com- 10 btcdcttwx lus noise environmem or a ^ fading envi . 
munication stations^When such a method breate down, say ronmem fof SDMA s ^ faave QDe ^ 
in low SINR or high fading situations, the method may not conventional channel, 
converge. This problem becomes more severe in communi- 
cation systems that have many users in the presence of high Another object is a adaptive smart antenna processing 
co-channel interference, that is, high interference from the 15 method and a PP aratus that determines adaptive smart 
signals within the conventional channel from other sub- anten , na P racessin 8 wei g hts efficiently in a low signal-to- 
scriber units when receiving a signal from a particular interference P lus noise environment or a high fading envi- 
subscriber unit, such other subscriber units being from the T0 * mc * for SDMA svstems mat have a P luraUtv of s P atial 
same or from neighboring cells in the case of a cellular channels per conventional channel, 
system that includes several receiving communications 20 Another object of the present invention is a adaptive smart 
stations, each communicating with a set of subscriber units antenna processing method and apparatus that determines 
located within its cell. adaptive smart antenna processing weights to use in a 
In theory, adaptive smart antenna processing permits current burst of data * the w «ghte adapted for the current 
more than one communication link to exist in a single burst of data by bemg determined from data from the current 
"conventional" communication channel so long as the sub- 25 burst. 

scriber units that share the same conventional channel can be ™,™ 

spatially (or spaUo-temporally) resolved. A conventional DESCRIPTION OF THE DRAWINGS 

channel includes a frequency channel in a frequency divi- The present invention will be more fully understood from 

sion multiple access (FDMA) system, a time slot in a time the detailed preferred embodiments of the invention, which, 

division multiple access (TDMA) system (which usually 30 however, should not be taken to limit the invention to any 

also includes FDMA, so to be precise, the conventional specific embodiment but are for explanation and better 

channel is a time and frequency slot), and a code in a code understanding only. The embodiments in turn are explained 

division multiple access (CDMA) system. The conventional with the aid of the following figures- 
channel is then said to be divided into one or more "spatial" FI(J x fc a filial bIock of a millt i.antenna 

channels, and when more than one spatial channel exists per 35 transceiver tem ^ indude a feceive wei ht deter . 

conventional channel the mulUplexmg is called space dm- mmer accQrdi tQ & of tQe , 
sion multiple access (SDMA). SDMA herein is used to mean _ T „ , , . , , „ 

inclusion of adaptive smart antenna processing both with . ™' 15 * ■ m0I f g ™ m ° f * transceiver 

one and with more than one spatial channel per conventional mclu ^g a ^ processor that when running a set of 

channel. 40 instructions implements a receive weight determiner accord- 

„ . j, . . . ... , ing to aspects of the present invention; 

Rapidly converging methods such as decision directed „ r . r 

methods also break down in the presence of high co-channel FIG - 3 B a flow chart of one embodiment of the weight 

interference in SDMA systems that have more than one determining method of the present invention; 

spatial channel per conventional channel. FIG- 4 shows a block diagram for a tracking reference 

Therefore, there is a need in the art for a adaptive smart 45 si 6 nal generator and demodulator used in the preferred 

antenna processing method that determines adaptive smart embodiment of the invention; 

antenna processing weights efficiently in a low signal-to- FIGS. 5 A, 5B, and 5C show the performance of, 

interference plus noise environment or a high fading envi- respectively, a constant modulus method, a decision directed 

ronment for SDMA systems that have one spatial channel method, and an implementation of the hybrid method 

per conventional channel, and for SDMA systems that have 50 according to aspects of the invention; 
a plurality of spatial channels per conventional channel. FIG. 6 shows a block diagram of the multiport weight 

Thus there is a need in the art for weight determination determiner and spatial processor according to the preferred 

methods that perform well under low SINR and high fading embodiment of the invention; 

situations, and that converge rapidly, i.e., in a small number 55 FIG. 7 is a flow chart of the preferred embodiment of the 

of iterations. multiuser weight determining method of the present inven- 

Thus there is a need in the art for a method that combines tion; 
good convergence properties with rapid convergence prop- FIG. 8 shows the effect of timing offset on the perfor- 

erties - mance of the CM method; and 

Thus there is a need in the art for a "blind" method (i.e., 6 o FIG. 9 shows a block diagram of an apparatus imple- 

one not using training data) that combines good convergence menting an aspect of the invention, 
properties (convergence when the SINR is low) with rapid 

convergence. DETAILED DESCRIPTION OF THE 



PREFERRED EMBODIMENT 
65 Base Station Architecture 

One object of the present invention is a weight determin- The preferred embodiment method and apparatus is 



SUMMARY 
esent inventio 

ing method that combines the advantages of methods that implemented in a communication receiver, in particular, in 
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a PHS-based antenna -array communications station embodiments described herein to incorporate less than full 

(transceiver) such as that shown in FIG. 1, with m antenna rate communication would be clear to those of ordinary skill 

elements in the antenna array. In the particular embodiment, in the art. 

m=4. While systems similar to that shown in FIG. 1 may be In the PHS used in the preferred embodiment, a burst is 

prior art, a system such as that of FIG. 1 with elements 5 defined as the finite duration RF signal that is transmitted or 

programmed or hard wired to carry aspects of the present received over the air during a single timeslot. A group is 

invention is not prior art. Also, this invention is in no way defined as one set of 4 TX and 4 RX timeslots. A group 

restricted to using the PHS air interface or to TDMA always begins with the first TX timeslot, and its time 

systems, but to any communications receiver that includes duration is 8x0.625=5 msec. 

adaptive smart antenna processing means. In FIG. 1, a to The PHS system uses jt/4 differential quaternary (or 
transmit/receive ("TR") switch 107 is connected between an quadrature) phase shift keying (n/4 DQPSK) modulation for 
m-antenna array 103 and both transmit electronics 113 the baseband signal. The baud rate is 192 kbaud. That is, 
(including one or more transmit signal processors 119 and m theriraitJ 192,0uTTsymbols per second, 
transmitters 120), and receive electronics 121 (including m FIG. 2 shows a more detailed but still simplified block 
receivers 122 and one or more receive signal processors 15 diagram of a PHS base station which includes adaptive 
123), for selectively connecting one or more elements of smart antenna processing and on which an embodiment of 
antenna array 103 to the transmit electronics 113 when in the the invention is implemented. Again, while systems with 
transmit mode and to receive electronics 121 when in the architectures similar to that shown in FIG. 2 may be prior 
receive mode. Two possible implementations of switch 107 art, a system such a that of FIG. 2 with elements pro- 
are as a frequency duplexer in a frequency division duplex 20 grammed or hard wired to carry aspects of the present 
(FDD) system, and as a time switch in a time division duplex invention is not prior art. In FIG. 2, a plurality of m antennas 
(TDD) system. The PHS preferred embodiment of this 103 is used, where m=4. More or fewer antenna elements 
invention uses TDD. The transmitters 120 and receivers 122 may be used. The outputs of the antennas are connected to 
may be implemented using analog electronics, digital a duplexer switch 107, which in this TDD system is a time 
electronics, or a combination of the two. The preferred 25 switch. When receiving, the antenna outputs are connected 
embodiment receivers 122 generate digitized signals that are via switch 107 to a receiver 205, and are mixed down in 
fed to signal processor or processors 123. Signal processors analog by RF receiver modules 205 from the carrier fre- 
119 and 123 may be static (always the same), dynamic quency (around 1.9 GHz) to an intermediate frequency 
(changing depending on desired directivity), or smart ("IF"). This signal then is digitized (sampled) by analog to 
(changing depending on received signals), and in the pre- 30 digital converters ("ADCs") 209. This is then downcon- 
ferred embodiments are adaptive. Signal processors 119 and verted digitally by digital downconverter 213 to produce a 
123 may be the same one or more DSP devices with different four-times oversampled complex valued (in phase I and 
programming for the reception and transmission, or different quadrature Q) sampled signal. Thus, elements 205, 209 and 
DSP devices, or different devices for some functions, the 213 correspond to receiver 122 of FIG. 1. For each of the m 
same for others. 35 receive timeslots, the m downconverted outputs from the m 

Note that while FIG. 1 shows a transceiver in which the antennas are fed to a digital signal processor (DSP) device 
same antenna elements are used for reception and 217 (hereinafter "timeslot processor") for further process- 
transmission, it would be cleat that one also can have ing. In the preferred embodiment, commercial DSP devices 
separate antennas for receiving and transmitting, and that are used as timeslot processors, one per receive timeslot. 
either only receiving or only transmitting or both receiving 40 The timeslot processors 217 perform several functions 
and transmitting may include adaptive smart antenna pro- including the following: received signal power monitoring, 
cessing, frequency offset estimation/correction and timing offset 

The Personal HandyPhone System (PHS), described for estimation/correction, smart antenna processing including 

example in the Association of Radio Industries and Busi- determining weights for each antenna element to determine 

nesses (ARIB, Japan) Preliminary Standard, Version 2, RCR 45 a signal from a particular remote user using a method 

STD-28 and variations as described in Technical Standards according to one aspect of the invention, and demodulation 

of the PHS Memorandum of Understanding Group (PHS of the determined signal. 

MoU — see http://www.phsmou.or.jp), is an 8 slot time divi- The output of the timeslot processor 217 is demodulated 

sion multiple access (TDMA) system with true time division data burst for each of the m (-4) receive timeslots. This data 

duplex (TOD). Thus, the 8 timeslots are divided into 4 50 is sent to host DSP processor 231 whose main function is to 

transmit (TX) timeslots and 4 receive (RX) timeslots. This control all elements of the system and interface with the 

implies that for any particular channel, the receive frequency higher level processing, which is the processing which deals 

is the same as the transmit frequency. It also implies with what signals are required for communications in all the 

reciprocity, i.e., the propagation path for both the downlink different control and service communication channels 

(from base station to users* remote terminals) and the uplink 55 defined in the PHS communication protocol. In the preferred 

(from users' remote terminals to base station) is identical, embodiment, host DSP 231 is also a commercial DSP 

assuming minimum motion of the subscriber unit between device. In addition, timeslot processors send the determined 

receive timeslots and transmit timeslots. The frequency band receive weights to host DSP 231. 

of the PHS system used in the preferred embodiment is RF controller 233 interfaces with the RF system, shown 

1895-1918.1 MHz. Each of the 8 timeslots is 625 micro- 60 as block 245 and also produces a number of timing signals 

seconds long. The PHS system has a dedicated frequency that are used by both the RF system and the modem. RF 

and timeslot for a control channel on which call initialization controller 233 receives its timing parameters and other 

takes place. Once a link is established, the call is handed to settings for each burst from the host DSP 231. 

a service channel for regular communications. Communi- Transmit controller/modulator 237, receives transmit data 

cation occurs in any channel at the rate of 32 kbits per 65 from host DSP 231. The transmit controller uses this data to 

second (kbps), called full rate. Less than full rate commu- produce analog IF outputs which are sent to the RF trans- 

nication is also possible, and the details of how to modify the miner (TX) modules 245, The specific operations transmit 
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controller/modulator 237 performs include converting data 

bits into a complex valued 7t/4 DQPSK modulated signal, ^ » (3) 

up-converting to a IF frequency, weighting by complex - 2j Ww(W-(*>- 

valued transmit weights obtained from host DSP 231, and ,=t 

converting the signals using digital to analog converters 5 

("DACs") to analog transmit waveforms which are sent to with spatio temporal processing, the convolution opera- 
transmit modules 245. lion of Eq. (2) is usually finite and on sampled data, 

Transmit modules 245 upconvert the signals to the trans- equivalent to combining the spatial processing with time 

mission frequency and amplify the signals. The amplified equalization using a time-domain equalizer with a finite 

transmission signal outputs are coupled to the m antennas io number of equalizer taps. That is, each of the w„{t) has a 

103 via dup lexer/time switch 107. finite number of values of t and equivalently, in the fre- 

Notatioo quency domain, each of the W ( (k) has a finite number of k 

The following notation is used. Let there be m antenna values. If the length of the convolving functions w„<t) is n, 

elements (m=4 in the preferred embodiment) and let ZjO), lhen ralner man determining a complex valued m-weight 

z^t), . . . , z m (t) be the complex valued responses (that is, 35 vector w r* one determines a complex valued m by n matrix 

with in-phase I and quadrature Q components) of the first, W r whose columns are the n values of w r (t). 

second, . . . , m'th antenna elements, respectively, after In the remainder of the description, whenever complex 

down-conversion, that is, in baseband, and after sampling valued receive weight vector w r or its elements are 

(four-times oversampling in the preferred embodiment). In mentioned, it will be understood that this may be for spatial 

the above notaTion7blit"riornecessarily for the invention, t is 20 processing or generalized to incorporate spatio-temporal 

discrete. These m time-sampled quantities can be repre- processing as described above to determine a weight matrix 

sented by a single m-vector z(t) with the i'th row of z(t) W r- Both spatial processing and spatio-temporal processing 

being z t (i). For each burst, a finite number, say N, of samples ^ thus referred to as adaptive smart antenna processing 

is collected, so that z^t), z^t), . . . , z ro (t) can each be herein, 

represented as a N-row vector and z(t) can be represented by 25 Determining Spatial Weights 

a m by N matrix Z. In much of the detailed description "Blind" methods for determining weights for adaptive 

hereinafter, such details of incorporating a finite number of smart antenna processing are methods that do not require 

samples are left out, and how to include these details would training data to be reconstructed. The methods of the present 

be clear to those of ordinary skill in the art. invention, like most blind methods, use some knowledge of 

Suppose several signals are sent to the base station from 30 lhe form of the originally transmitted signal and constrain 

several, say N„ remote users. In particular, suppose that a the 0Ut P ut si &* 1 t0 have one or more ^ si ^ x 

subscriber unit of interest transmits signal s(t). Adaptive properties. The property may be the amplitude 

smart antenna processing includes taking a particular com- characteristics, or some statistical characteristics such as the 

bination of the I values and of the Q values of the received entropy or cyclo-stationarity, the correct modulation 

signals z^O, z^t), . . . , zjt) in order to extract an estimate 35 scheme, or reconstructing an accurate replica. Such methods 

of the transmitted signal s(t). Such weights may be repre- are sometimes called "property restoral." 

sented by the receive weight vector for the this particular Methods With Good Convergence Properties 

subscriber unit, denoted by a complex valued weight vector Methods with good convergence properties include par- 

w r with i* element w„.. The estimate of the signal sent is lial P ro Perty restoral methods: those methods that reproduce 

then 40 one or more properties without attempting to reproduce an 

exact replica by determining the bit stream and reconstruct - 

(I) ing a signal. In this class are included methods that preserve 

3(r) = £<z,(0 = the amplitude (modulus), the entropy, and the spectral 

- i=I coherence (e.g., cyclo-stationarity) of the signal. 

45 The constant modulus (CM) method is a very simple and 

, . ... , . £ h • L effective technique applicable for signals modulated by 

where w „• is the complex conjugate of w„. and w r » is the schemes that resul , in constantamrjtui^^ These 

Henmlian transpose (that the transpose and complex include a „ forms of phSTTnrfquentyl^ulation, 

conjugate of receive weight vector w r In embodiments mcludi , he differenlial hase shift keying modulation of 

which include spatio-temporal processing each element in so the pHS tem ^ m , he hmd embodiment . ^ 

the recede we.ght vector ,s a funct.on of time, so that the lained be , me CM method ako fa ^ tQ Qon 

weight vector may be denoted as w^t), with ,th element conslant moduhls si ^ ^ CM me *£ d determines 

w„<t). The estimate of the signal may then be expressed as weights ^ a) restore the (constant 

modulus) property of a signal and/byproduce a signal that, 

3(0 = Y W (0*c (0 (2) 55 ^ transmitted by a remote user, would produce signals at the 

& " antenna elements of the array that are "close" to the signals 

actually receivedfAmplitude variations may be introduced 
by interference, fi^-fading, and by timing offsets in those 

where the operator " *" is the convolution operation. Spatio- nioclulationschemes wherein the constant modulus property 

temporal processing, for example, combines time equaliza- 60 depends on accurate timing offset correction,]including, for _ 

tion with spatial processing, and is particularly useful for example, the DQPSK modulation scheme Trf the preferred 

wideband signals. Forming the estimate of the signal using embodiment where the constant modulus property only 

spatio-temporal processing may equivalently be carried out holds at baud points. In the presence of co-channel 

in the frequency (Fourier transform) domain. Denoting the interferes, the CM method tends to pick the strongest 

Frequency domain representations of §(t),z ( <t), and w„<t) by 65 signal, whether this is the desired signal or a co-channel 

§(k), Z/k), and W t <k), respectively, where k is the discrete interferer. Even if the desired signal strength is greater than 

frequency value, any interferer's by only 0.5 dB, the CM will still correctly 
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pick up the strongest, i.e., the desired signal. That is, a CM 2. For the samples of interest, perform a copy signal and 

method has very good convergence properties. normalize: 

There are many variations of the CM method. These 

typically minimize a cost function of the general form: ( 5 ) 



■WO = 



= £)E|5^(0l p -in (4) 



(6) 



w r"R a ~ X r„ f (7) 

where R =ZZ" r « 



3. Compute the weight vector w r using a least squares 

where E(.) denotes a statistical expectation operation, and p procedure. That is, 

and q are positive integers, typically 1 or 2. It would be clear io 

to those in the art that, in practice, the statistical operation is N 

replaced with some form of sample averaging or accumu- w r = argmin^T (5,^(0- 

lation (e.g., by summation for a set of samples, which in the Wr *=i 
preferred embodiment is a subset of all the samples in a 

burst). Also, it would be clear and is within the scope of the 15 

present invention to add more terms in the cost function of where N fe the number of samples used in the calcu 
Eq. (4). For example, a term may be added to limit the lation * ^ solution of Eq. (6) is 
weight vector magnitude. See above-referenced U.S. appli- 
cation Ser. No. 08/729,390 for an example of a cost function 
(not a CM copst function) with such an added term. Th 20 
signal s w y(t) is the normalized copy signal (called the 
"reference signal") used in the cost function. That is, the 
reference signal for determining the weights is the weighted r - V co- 
sum of received antenna signals which is then normalized. "it ™* 
Weight determining is determining the set of weights that 25 
minimizes the cost function in Eq. (4). 

The CM method also is applicable to non constant modu- z(t)s' nf (t); and N is the number of samples used; and 

lus signals. See, for example, J. Lundell and B. Widrow: 4. Repeat steps 2 and 3 until convergence is reached, 

"application of constant modulus adaptive beamformer to Note that in the calculation of step 3, in practice, overall 

constant and non-constant modulus signals," Proceedings, 30 scale factors are unimportant. All the scale factors for the 

1988 Asilomar Conference on Signals, Systems and Com- weights preferably are applied in combination as a gain in 

puters (AGSSC-1988), pp. 432-436, 1988. Lundell and the system. 

Widrow use a cost function like that of Eq. (4) with p=q=2 Note that the CM method can be extended to spatio- 
(this is called a 2-2 CM method) and showe that any constant temporal processing. One known method uses a 2-2 CM 
and non-constant modulus signal can be recovered using 35 method and shows that under certain assumptions usually 
such a 2-2 CM method, as long as the ratio of the fourth met in practice, the CM method for spatio-temporal weight 
moment to the square of the second moment (the ratio called determination (i.e., for weight matrix determination) is sure 
the kurtosis) is below 2. For example, an M-quadrature to converge. See C. B. Papadias and A. Paulraj, a space-time 
amplitude modulated signal (M-QAM), is known to have a constant modulus algorithm for SDMA systems," 
kurtosis of approximately 1.4-1. 2/(M-l), hence the kurtosis 40 Proceedings, IEEE 46th Vehicular Technology Conference, 
of any QAM signal is always smaller than 1 .4. A CM method pp. 86-90, 1996. The Papadias et al method, however, is not 
therefore is applicable to such a signal. block data based. However, the spatial weight determining 
At least one iteration of a particular method with good method can easily be modified for spatio-temporal process- 
convergence properties, the CM method, is used in the ing according to a weight matrix by re-expressing the 
preferred embodiment, and that implementation of the CM 45 problem in terms of matrices and vectors of different sizes, 
method uses values of 1 for p and 2 for q in Eq. (4). When As throughout this description, let m be the number of 
this invention is applied to a non constant modulus signal, antenna elements, and N the number of samples. Let n the 
and the CM method is used, other values for p and q, e.g., number of time equalizer taps per antenna element. Each 
p=q=2, may be used. The preferred implementation is also row vector of N samples of the (m by N) received signal 
a block based method. That is, a block of the antenna 50 matrix Z can rewritten as n rows of shifted versions of the 
received signals is weighted and the weights are determined first row to produce a received signal matrix Z of size (mn 
using this block of data. The block is a subset of the samples by N), which when pre-multipiied by the Hermitian trans- 
in a burst. In particular, 75 samples of the 120 PHS burst pose of a weight vector of size (mn by l)produces an 
symbols preferably are used, where the 75 symbols are in the estimated received signal row vector of N samples. The 
payload which is in the middle of the PHS burst. Using data 55 spatio-temporal problem has thus been re-expressed as a 
from the payload advantageously ensures that the data used weight vector determining problem. For the CM method, in 
for weight calculation for any one remote user is not the Eq. (7), for example, the weight vector is a "long" weight 
same as that for another subscriber unit. There are a maxi- vector of size (mn by 1), is a matrix of size (mn by mn), 
mum of 88 such payload samples in a PHS burst. and r^ is a long vector of size (mn by 1). Rearranging terms 

With p-1 and q-2, the method is called the least squares 60 provides the required (m by n) weight matrix, 

constant modulus method, and includes the following steps. Since in the preferred embodiment, in order for the CM 

1. Initialize the weight vector. For example, use property to hold, the sampled data needs to be approximately 

w r,*mxw/~[l 0 0 . . . 0J where x* denotes the transpose on-baud, in carrying out step 2, timing offset estimation and 

of x. In an improved embodiment, the largest eigen- correction are performed, which in this case may include 

vector of R^-ZZ" corresponding to the largest singular 65 decimation in time and interpolation because the samples of 

value of Z is used. In yet another embodiment, the the received signals from each antenna element of antenna 

weight vector from the previous burst is used; array 103 are oversampled and may include some timing 
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offset. The variable t in Eqs. (5) and (7) therefore represents 
approximately on-baud times for the samples. As would be 
clear, and as is also described in above-referenced U.S. 
patent application Sen No. 09/153,110, the timing offset 
estimation/correction (which may include decimation/ 
interpolation) may be performed on the m signals prior to the 
signal copy operation, or on the signal after the signal copy 
operation. 

Simulations were performed to determine the level of 
accuracy of the timing offset/baud point estimation. FIG. 8 
shows the results. In the simulation, the output SINR 
obtained using CM weights was calculated when the signal 
was sampled not exactly at the baud point, but offset from 
the ideal baud point by some timing offset which varied from 
- V£ to +Vi baud in steps of Y& of a baud. The results show that 
even for a signal offset by ±Vs of a baud, the output SINR 
degrades only by 0.6 dB. Although this number is specific to 
the test case shown, the conclusion is that the accuracy of the 
timing offset correction for the CM method need not be high. 
A simple method therefore may be used for the offset 
correction/decimation/interpolation to produce approxi- 
mately baud-aligned samples. 

Also note that timing offset correction (including 
decimation/interpolation) of an oversampled signal may not 
be necessary for all modulation schemes that have the 
constant modulus property. For example, the Digital Euro- 
pean Cordless Telecommunications (DECT) standard and 
the Global System for Mobile Communications (GSM) 
standard use Gaussian Minimum -Shift Keyed (GMSK) 
signals, which always have a constant modulus, so timing 
offset correction is not needed for the CM method in those 
cases. 

Implementation of a least squares CM weight determina- 
tion method is fairly straightforward. Since no demodulation 
occurs, frequency offset estimation and correction, and 
demodulation, are not required. In one embodiment of the 
present invention, frequency offset correction is carried out 
when implementing the CM method, even though this is not 
required. An additional feature of simple property restoral 
methods such as constant modulus methods is that conver- 
gence occurs even under very low signal-to-interference- 
plus-noise (SINR) values. 

The main disadvantage of methods with good conver- 
gence properties such as simple property restoral methods is 
that they may take many iterations to converge. In a typical 
system, processing power, such as DSP processing power, is 
extremely limited, and thus convergence using the CM 
method may not occur within a specific amount of time, for 
example in order to use the weight vector within the current 
burst. 

Rapidly Converging Methods 

Rapidly converging methods such as decision directed 
methods, unlike partial property restoral methods, converge 
very fast. With a decision directed method, the property 
restored is a complete replica of the originally transmitted 
signal that has the correct modulation scheme. That is, a 
signal copy operation, such as Eq. (1), estimates a received 
signal, and that signal is demodulated and a reference signal 
with the correct bit stream is constructed. To work well, one 
needs to correct for any frequency and timing offsets when 
constructing the reference signal. The correct weights are 
those that produce a reference signal that is close to the 
transmitted signal. The scheme may involve one or more 
iterations to achieve the "best" weights. While the timing 
offset correction (including any decimation) and the fre- 
quency offset correction are shown below as occurring after 
the signal copy operation, clearly one or more of these can 
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be carried out prior to the signal copy operation. See 
above-referenced co-owned U.S. patent applications Ser. 
Nos. 08/729,390 and 09/153,110 for examples of these 
operations occurring before and after signal copy, and for a 
detailed description of a decision directed method. When a 
least squares criterion is used, the general method includes 
the following steps: 

1. Initialize the weight vector. For example, use 
W r>una/-{1. 0 0 ... 0]' where x* denotes the transpose 
of x. In an improved embodiment, the singular vector 
of R^TZ" corresponding to the largest singular value 
is used. In yet another embodiment, the weight vector 
from the previous burst is used. As will be described 
below, one aspect of the present invention includes 
using a decision directed method after a partial property 
method is used. In such a case, when implemented in 
any of the embodiments of the present invention, the 
last obtained weight vector (i.e., using a partial prop- 
erty restoral method) is used; 



55 



60 



65 



2.; Perform a signal copy 



J<f)-HV"z</) f 



(8) 



followed by decimation/interpolation if the samples are 
originally oversampled (in an alternative, the 
decimation/interpolation may occur prior to the copy 
\signal operation); 

3. Estimate llrnrng-^fTd frequency offset to produce a 
signal that has the correct timing and frequency offsets; 

4. Determine a reference signal s rc J(i) by making symbol 
decisions (i.e., demodulating), such that s re ^(t) has the 
correct bit stream and the same modulation scheme, 
and the same timing and frequency offsets as the signal 
transmitted to the receiver from the particular user; 

5. Computing the weight vector by least squares minimi- 
zation of over w r . That is, 



w r = arg min £ \s^(t) - y^!z{t)\ 2 



which has the solution 



where R -ZZ and 



(9) 



(10) 



8 



and 



6. Repeat steps 2, 3, 4, and 5 until convergence is reached. 

Note that steps 2, 3 and 4 require that the signal be 
corrected for frequency and timing offsets so that the correct 
demodulation decision is made in step 4, while step 5 
typically requires that the correct frequency and timing 
offsets are re -introduced so that the reference signal and the 
copy signal in the cost function have the same timing and 
frequency offsets. Note also that the minimization of Eq. (9) 
also may include other terms, such as a weighted norm of the 
weight vector term to place a constraint on the norm of the 
weight vector, as described above for the CM case and as 
described in above-referenced U.S. application Ser. No. 
08/729,390. See also above-referenced co-owned U.S. 
patent applications Ser. Nos. 08/729,390 and 09/153,110 for 
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a detailed description of how to determine a reference signal then determined using averages over the number of samples 

(step (4)). This also is described below with the aid of FIG. for the expectation operation. 

4. 

Note that the decision-directed method can easily be j * 00 

extended to determining the weight matrix for spatio- 5 R = — ^j^tO+C^tO.and 

temporal processing, for example, by re-arranging terms as r=l 

described hereinabove for the CM method, and by other ^ N (12) 

methods as would be clear to those in the art. This invention, ff 1 = — ^ (/ 2 W + Q 2 ^)) 2 . 
therefore, also covers methods for determining weight vec- 

tors and weight matrices for spatio-temporal processing. 10 

Decision directed methods thus reproduce an accurate Nole lhat once me instantaneous power R^oPO+Q^t) is 

replica of the signal supposedly transmitted to the receiver, determined, determining the squared power R 4 (t>[R 2 (t)] 2 

whereas partial property restoral methods reproduce one or requires only a single additional multiplication per sample, 

more simple properties, such as the correct amplitude. and the estimated signal-to-interference-plus-noise-ratio is 

Decision directed systems perform extremely well, and 15 determined, preferably with at most one square root 

converge in very few iterations in a reasonably high SINR operation, using 
environment. However, these methods are sensitive to the 
initial conditions, and may not even converge when the 
initial SINR is low. This is common in high mobility cellular 



— ~ (12) 
2-& 



systems and other systems that exhibit fading. 20 ^ (K*)' 

Note that an iteration of the CM method is typically sinr- 
computationally cheaper than an iteration of a decision t _ 



directed method, not requiring frequency offset correction v (fif 

nor demodulation. 

The Preferred Method: Single User 25 _ a-Ja ^ r 4 

One aspect of the invention is a weight determining ~ 1 - a * w m ^2 
method which includes carrying out a number N 2 of itera- 
tions of an iterative weight determining method with good 

convergence properties, such as a partial property restoral Both the ratio 

method, preferably the CM method, combined with and 30 
followed by a second number N 2 of iterations of a rapidly 

converging method, such as the decision directed method to ~~i 

achieve the advantage of good convergence properties with ^ R ^ 
rapid convergence. The N a CM iterations bring the starting 

condition for the N 2 iterations of the decision directed 35 and the quantity A are sometimes called the kurtosis. This 

method to the region where rapid convergence of the deci- preferred method of signal quality estimation is insensitive 

sion directed method is most certain. The preferred embodi- to frequency offset, so is a particularly attractive method for 

ment uses one iteration of the decision directed method use with the CM method which also is insensitive to 

(N 2 =l) while an alternate implementation uses two frequency offsets. 

iterations(N 2 =2). The method can be re-stated as carrying 40 Other methods of determining the quality of the post copy 

out iterations of an iterative weight determining method with operation signal also may be used in alternate implementa- 

good convergence properties until a switching criterion is tions of the invention. 

met, then, starting with the weights obtained with the The method for weight determination for a single user is 

method with good convergence properties, carrying out a illustrated in the flow chart of FIG. 3. An initial weight 

number of iterations of a rapidly converging method. In one 45 vector is formed 303. This might be [1 0 0 . . . 0]', or in an 

case, the switching criterion is an explicitly defined number improved embodiment, the singular vector of R^-ZZ 7 * 

N x of iterations. In another, the preferred embodiment, Nj is corresponding to the largest singular value is used. In yet 

not explicitly specified. Rather, the switching criterion is a another embodiment, the weight vector from the previous 

SINR threshold for the copy signal, and the switchover to burst is used. A copy operation 305 is now performed 

the decision directed method occurs when the SINR estimate 50 according to Eq. (1), but in the preferred embodiment, using 

is equal to or exceeds the threshold. In this way, a sufficient only the middle part of a burst, preferably only 75 symbols 

number Nj of CM iterations is used in order to achieve a (300 samples) from the payload part in the middle of the 

SINR which is sufficient to ensure convergence of the burst. The output is corrected for timing offset 307. Any 

decision directed method in only N 2 further iterations. timing offset correction method may be used. As discussed 

Many methods of determining SINR estimates may be 55 hereinabove, the timing offset correction need not be very 
used. In the preferred embodiment, the method used is as accurate. The preferred method is as described in above- 
described in U.S. patent application Ser. No. 09/020,049 to referenced U.S. patent application Ser. No. 09/153,110. The 
Yun entitled POWER CONTROL WITH SIGNAL QUAL- copy and timing offset correction operation may be com- 
ITY ESTIMATION FOR SMART ANTENNA COMMU- bined. Inherent in the timing offeet correction operation, but 
NICXTION SYSTEMS (filed Feb. 6, 1998). The implemen- 60 not explicitly shown in FIG. 3, is any necessary decimation 
tation of the signal quality estimation method is now and interpolation, so that after step 307, the data includes 75 
described. complex valued (I and Q) samples approximately at the baud 

Denote by N the number of samples of a burst to use for points of 75 symbols from the middle of the current burst, 

the estimate. The sampled modulus information is first The SINR of the copied signal is estimated 309 preferably 

extracted by forming the sum of the squares of the in phase 65 using the kurtosis as described hereinabove. In step 311, it 

and quadrature signals (the real and imaginary parts of is determined if the SINR has exceeded a threshold SNR. If 

signal s(t). The mean power and mean squared power are not, an iteration of the constant modulus method is per- 
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formed in step 313 using the minimum squared cost function 
criterion as described by Eqs. (6) and (7) above. Then the 
method returns to the copy operation of step 305 for another 
iteration. If on the other hand, the SINR threshold is deter- 
mined in step 311 to have been exceeded, N 2 iterations of the 
decision directed method are performed in steps 315 and 
317, including frequency offset correction 315. Any fre- 
quency offset correction method may be used, and the 
preferred method is as described in above-referenced U.S. 
patent application Ser. No. 09/153,110. Similarly, for the 
decision directed adaptation, including generating the ref- 
erence signal, any method may be used, and the preferred 
embodiment uses the method described in above-referenced 
U.S. patent application Ser. No. 09/153,110. When the 
weight is determined, in the preferred embodiment, only a 
subset of the samples in each burst is used. Therefore the 
finally determined weight vector is now used in a copy 
operation and demodulation step 318 on the whole burst. In 
this embodiment, step 318 includes timing and frequency 
offsets determination and correction, demodulation, prefer- 
ably using the architecture described below with the aid of 
FIG. 4. The output of the decision directed adaptation is a 
signal 319. 

The preferred embodiment for the reference signal gen- 
eration which is part of decision directed adaptation step 317 
(using part of the burst data), and is used for demodulating 
all the burst data in step 318, preferably uses a reference 
signal generation architecture and method that includes a 
tracking mechanism, preferably sample-to-sample, that 
forms the phase of the reference signal at a sample point by 
relaxing the phase of a signal ideally advanced from a 
previous reference signal sample, towards the phase of the 
copy signal at the same sample point, the copy signal formed 
from the received antenna signals. The reference signal is 
constructed at each sample point, by constructing an ideal 
signal sample from the copy signal at the same sample point, 
the ideal signal sample having a phase determined from the 
copy signal at the sample point, with the phase of the ideal 
signal sample at an initial symbol point set to be an initial 
ideal signal phase, and relaxing the phase of the ideal signal 
sample towards the copy signal sample phase to produce the 
phase of the reference signal. The phase of the ideal signal ' 
is determined from the phase of the reference signal at the 
previous sample point for which the phase is determined, 
and from a decision based on the copy signal. In one 
implementation, the reference signal is determined in the 
forward time direction, and in another implementation, the 
reference signal samples are determined in the backwards 
time direction. In one version, the step of relaxing the phase 
of the ideal signal sample towards the phase of the copy 
signal b^n) corresponds to adding a filtered version of the 
difference between the copy signal phase and ideal signal 
phase. In another version, the step of relaxing the phase of 
the ideal signal sample towards the phase of the copy signal 
corresponds to forming the reference signal sample by 
adding to the ideal signal sample a filtered version of the 
difference between the copy signal and ideal signal. 

This is now described in detail with the aid of FIG. 4 using 
the ji/4 DQPSK PHS signal as an example. Modifying for 
other modulation schemes would be clear to one in the art. 
Phase detector unit 403 detects the phase difference 405 
between the copy signal 401 (corrected for the timing and 
frequency offsets) and the previous reference signal 417. 
The phase difference signal 405 is fed to a slicer 407 to 
generate the decision phase difference 419. The correct 
phase difference for jt/4 DQPSK is (2i-l)*/4, i-1, 2, 3, or 
4, and is the phase difference between the previous reference 
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signal sample and the ideal signal. This is subtracted in block 
409 from the actual phase difference 405 to generate error 
signal 411. This error signal is filtered in filter 413 to 
generate filtered error signal 415. It is this filtered error 
signal that is used to adjust the phase difference 419 closer 
to the actual phase difference 405. The corrected phase 
difference 421 is then used in a frequency synthesizer/phase 
accumulator 423 to generate the reference signal 429. It is 
the previous sample value 417 of the reference signal 429 
that is used by phase detector 403, so a unit time delay 425 
is shown between these signals. The symbols of the signal 
430 (signal 319 after the N 2 iterations) are determined by 
block 427. Mathematically, if b^t) denotes the reference 
signal complex sample values at the baud point t, and L 
denotes the phase, the input to phase accumulator 423, 
^(t)-^(t-l), is 

Wtcr{Z.d We(I Xn)-decidc{£d i(/M Xn)}}+dccide{£d WM Xn)}, 

where decide {Z.d w<ra Xn)} is the output of slicer 407 and 
equals (2i-l)jtf4, 2, 3 or 4 for k/4 DQPSK. In this, the 
"ideal" complex valued sample points b idea £l) are defined as 
follows: 

WOW>*(oH>ro), 

where b(t) are the samples of input signal 401, and ^.d^/i) 
is the phase difference between the present input sample and 
the previous reference signal sample: 

^WO-^0-^('-i)=4WV(M)l 

where * denotes the complex conjugate: The "ideal" signal 
is the reference signal with the phase advanced by an ideal 
amount, that amount depending on the decision made 
according to £d lrfeo/ (0- 

^«X0=^(f-l)+(2*-l)*/4> W, 2, 3, or 4. 

To obtain the reference signal, the phase of b^XO k now 
relaxed towards the phase of b(t) by filtering the quantity 
[Z.b(t)-Z.b W(ra X0]» tne phase error between b(t) and b idee Jl), 
and adding the filtered quantity to the phase of b^XO- An 
alternate embodiment filters the quantity (b(t)-b w<fa X0) 
rather than the phase error. The filter preferably is a constant 
of proportionality. Higher order filters may be sued. 
Mathematically, in one embodiment, 

and in another embodiment, the architecture of FIG. 4 may 
be modified slightly to use 

In the preferred embodiment, the method of the flowchart 
of FIG. 3, including the tracking reference signal generator, 
is implemented as a set of instructions for the timeslot 
processor 217 which is a signal processor (DSP) device. 

A simulation of the method of FIG. 3 was carried out for 
the system described, but with the initial weight vector the 
eigenvector of R„oZZ" corresponding to the largest eigen- 
value. The simulation was carried out for the PHS base 
station of FIG. 2 with four antenna elements. The input 
signals at each of the antennas had a signal-to-noise ratio 
(SNR) of 11.9 dB. The input carrier to interference ratio 
(CIR) was 1.1 dB, corresponding to an initial copy signal 
SINR of 0.8 dB. While a regular PHS burst has 120 symbols, 
only the center 75 symbols were used for the weight 
calculation. All calculations were carried out offline using 
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the MATLAB environment (The Mathworks, Inc., Natick, 
Mass.)- The results are shown in FIGS. 5A, 5B, and 5C, 
which compare the convergence characteristics of the deci- 
sion directed method, the least squares CM method, and the 
combined method of the present invention under low SINR 5 
conditions (in this case, a post-copy SINR of about 0.8 dB 
after the first copy operation with the initial weight vector). 
The output SINR (dB) after each iteration is measured and 
plotted using the SINR estimation method. The first SINR 
value shown is when the SINR is estimated after the copy 10 
operation with the initial weight vector (first singular vector 
of RJ. This is the same for all three methods. As seen on 
FIG. 5 A, the decision directed method does not converge 
even after 10 iterations. FIG. 5B shows that the CM method 
converges slowly, with the output (estimated) SINR steadily 

^ re a^fh ^^l^^ in°a ^aK^s^ m 0 P^ ma ^ SINR is 18 
converge to this optimal value. FIG. 5C shows the method 
of the present invention operating with a switchover output 
SINR threshold of 7.5 dB. Notice that the results start off 
identical to that of FIG. 5B but then diverge after the 
switchover to the decision directed method (the results of 
FIG. 5B are shown dashed after the switchover). After the 
decision directed method starts, the method reaches the 
optimal SINR after only 2 iterations of the decision directed 25 
method, and gets very close to the optimal even after only 
one iteration of the decision directed method. Overall, the 
method of the flowchart of FIG. 3 converges within 5 
iterations as opposed to more than 10 iterations if only the 
CM method is used. 
Multi-port Architecture 

When there are several, say N^, subscriber units, both 
in-cell and out-of-cell, that are in th e same conventional 
^ channel (i.e., co-channel users), ^he preferred embodiment 
of this invention uses a multi-port architecture, each "port" 
separately forming the copy s ignal and tracking a single one 
ol^He^NJljubjiuibLi-uMils kuy^port's subscriber unit thus 
being a co-channel interferer to the other subscriber units 
and their corresponding ports. Only co-channel users that 
*^rTave signal components received at the antenna elements 40 
that are above some noise floor are so tracked. The number 
of such users may be estimated. Given any burst (matrix Z), 
the eigenvalues of R^ZZ" may be examined and an order 
estimation may be performed. Any order estimation method 
may be used. For example, the Rissanen minimum descrip- 
tive length (MDL) criterion, or the Akaike information 
theoretic criterion are well known. For a survey of tech- 
niques for determining the number of significant co-channel 
users, see Section 3.8 of Rias Muhamed and T. S. Rappaport: 
"Direction of arrival estimation using antenna arrays," Tech- 
nical Report MPRG-TR-96-03, Mobile and Portable Radio 
Research Group, Bradley Department of Electrical 
engineering, Virginia Polytechnic Institute, January 1996 
(also Rias Muhamed: "Direction of arrival estimation using 
antenna arrays," Masters Thesis, Bradley Department of 55 
Electrical Engineering, Virginia Polytechnic Institute and 
State University, Blacksburg, Va. 24061, USA). The pre- 
ferred embodiment uses the minimum descriptive length 
criterion. 

While the preferred embodiment includes estimating and 60 
then tracking all the "significant" co-channel users, in 
another embodiment, a "good" wireless design environment 
is assumed, that is, co-channels users not communicating 
with the same base station are assumed far away, so that the 
. only co-channel users of significance are those users that 65 
share the same conventional channel and communicate with 
the base station. That is, the subscriber units that are the 
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different spatial channels in the conventional channel. In 
such a case, N, is known. 

Considering, for example, when there are two known (by 
estimation or by knowledge) co-channel users (i.e., N,=2) in 
the environment. When tracking one of the subscriber units, 
the other subscriber unit is an interferer. Therefore in this 
architecture, the desired signal and the significant interferes 
communicating in the same conventional communication 
channel are simultaneously tracked. In a fading environment 
as encountered, for example when the subscriber units move 
rapidly, the carrier-to-interference-ratio (CIR) may be very 
low and the instantaneous CIRs can fluctuate over a wide 
range. At any point in time, at any given burst, the desired 
signal therefore may be weaker than any of its interferers 
and the desired signal's port may become interferer-locked. 
That is, it may start tracking an interferer rather than a 
desired remote user. 

FIG. 6 shows the block diagram of the preferred embodi- 
ment multi-port adaptive smart antenna processing appara- 
tus. In each port, t he oversampled outputs 6 05 of receivers 
122 from the antenna elements 1U3 are combined in a signal 
copy operation 607, initially using an initial weight vector 
631 -i, i=l, . . . , Ny for the first, .... N 5 th port, respectively, 
these initial weights provided by a weight initializer 621. 
The resulting copy signal is timing offset corrected by timing 
offset corrector unit 609 which also decimates/interpolates 
to produce a set of approximately baud-augnecTsamples (for 
the CM method iterations) or substantially baud-aligned 
samples (for the decision directed method iterations)) baud- 
aligned samples. The baud-aligned samples are fed into an 
SINR estimator 613 the output of which is fed into a weight 
calculator and demodulator 615 which uses the baud-point 
aligned samples and/or the antenna signals 605 to determine 
a reference signal and a set of weights according to the 
inventive method described herein. Since at least one itera- 
tion of a decision directed method is used in weight calcu- 
lator and demodulator 615, the output is a demodulated 
signal 617. In this way, N s demodulated signals for N r 
subscriber units are determined. Having multiple ports 
enables tracking both any desired subscriber unit-sent 
signal, and any c p-channel interferer/ The adaptation method 
further described below nas the ability to switch between 
any desired user and th e interferer in a fading environment. 
So by using N s ports, the N s users are tracked at the same 
time, and if any user's signals jumps from one port to the 
other, which might occur in a fading environment, the output 
of the ports is sorted in user sorter 623 to separate the desired 
user from a ny interferers . and correctly output N, demodu- 
lated signals '62^ 

Other multiport architectures are known but not for use 
with the adaptive methods described herein. See, for 
example, B. G. Agee: "Blind separation and capture of 
communication signals using a multitarget constant modulus 
beamformer," 1989 IEEE Military Communications Confer- 
ence ("MILCOM 89"), vol. 2, pp. 340-346, New York: 
IEEE, 1989, for a multiport architecture for a constant 
modulus method. The Agee method differs in several 
respects from the method described herein, including for 
example, 1) in the weight initialization method; 2) in what 
computations are carried out in each of the ports. Agee's 
method jointly orthogonalizes all the weight vectors at each 
stage of the iteration, which is computationally expensive, 
while the preferred embodiment of the present invention, 
after jointly initializing each of the ports, allows each port to 
adapt independently; 3) the method of determining weight 
vectors is different. Note that since it is anticipated that in the 
future, computational power will become more readily 
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available, in an alternate embodiment, each port's weight 
vector is made orthogonal. 

The Preferred Weight Determining Method: Multiple User 
FIG. 7 describes by means of a flow chart preferred 
method for determining the weights and the output signals 
625. To start with, the initial copy is performed using the 
eigenvectors of the R a matrix. The port labeled "#1" is 
initialized in step 703 with the eigenvector 631-1 corre- 
sponding to the largest singular value, the second port with 
the next eigenvector 631-2, . . . , and the N f th port with the 
Njth eigenvector 631-N,. These eigenvectors are guaranteed 
to be linearly independent, and also, are in general a pre- 
ferred value to start with. Alternatively, any substantially 
independent initial weight vectors may alternatively be used. 
For example, in an alternate embodiment, port #1 is initial- 
ized with the vector [10 0... 0] ! , port #2 with the vector 
[0 1 0 0 ... 0J, and so forth. In each port, after initialization, 
the method proceeds in the manner in each port as for the 
single user case of the flow chart of FIG. 3. That is, step 305 
is the copy operation performed first with the initial value. 
The resulting signal is in step 307 timing offset corrected 
(which includes decimation/interpolation if initially 
oversampled) to produce substantially baud aligned samples 
which are fed into a signal quality estimator which in step 
309 estimates the SINR substantially on the baud points. A 
decision is made in step 311 whether to choose the CM 
method or the decision directed method for the weight 
adaptation. If the SINR is below a predefined SINR 
threshold, an optimization based on a partial property resto- 
ral method (preferably the CM method) is carried out in step 
313 and the method returns to step 305 for another iteration, 
starting with the most recently determined weight vector for 
this port. If the SINR is above the threshold, then in step 315 
frequency offset correction is carried out, and a single 
decision directed adaptation iteration is carried out in step 
317. Note that if the timing offset correction was only 
approximate in step 307, a more accurate correction would 
be needed for the decision directed method, and such a 
modification would be clear to one of ordinary skill in the 
art. Only one decision directed iteration is carried out in the 
preferred embodiment. Alternatively, more than one deci- 
sion directed iteration is carried out. When the weight is 
determined, since in the preferred embodiment, only a 
subset of the samples in each burst is used, the finally 
determined weight vectors for each port are used in a copy 
operation and demodulation step on the whole burst. This 
copy operation in the preferred embodiment includes timing 
and frequency offsets determination and correction, and 
demodulation, the demodulation preferably using the archi- 
tecture shown in FIG. 4. The result for each port is a 
demodulated signal 617. 

Note that one feature of the invention in both the single- 
user and multi-user case is using the weight obtained using 
the current burst's data for determining the signal for the 
current burst. When the subscriber units move around, and 
in other fading and low SINR environments, using a weight 
vector from the previous burst may not produce good results. 

The final step is the sorting of these outputs to determine 
if any of the output ports have become interferer-locked. 
PHS bursts, for example traffic channel bursts, include fields 
for the payload, a unique word (UW) that is known to all 
subscriber units, and a error detecting cyclic redundancy 
check (CRC) field. Other protocols include somewhat dif- 
ferenrtelds iharmay be used to determine if a particular 
message is from a particular subscriber unit or for a par- 
ticular base station. To determine interference locking, in the 
preferred embodiment, the idea is to distinguish a desired 
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subscriber unit sending a waveform valid for the system 
from an interfering subscriber unit also sending a waveform 
which is valid for the system. That is, means exist for 
defining a "valid" subscriber unit waveform, for example, 
the waveform having some required data and modulation 
format and scrambled with a particular key for the sub- 
scriber unit. An interfering unit similarly would include 
some means for defining its own validity, for example the 
waveform having some required data and modulation format 
and scrambled with a particular key different from that for 
the subscriber unit. In the preferred PHS implementation, 
one method for detecting interferer-locking includes moni- 
toring both the unique word (UW) and the CRC. It is 
specified that the data bits in each burst be scrambled with 
a bit pattern that is generated using the lower 9 bits of the 
cell station identification code (CSID). This 9-bit word used 
to cipher both the burst payload and its associated CRC is 
called the scrambling key. When designing a communica- 
tions system, for example a cellular system, it is advisable to 
make sure that neighboring communications stations (base 
stations) each have a different scrambling key. Note that in 
the PHS specification, a base station or communications 
station is called a cell station. 

Knowing that, the preferred embodiment basic interferer- 
locking detection method involves the following steps: 

For a particular port, that is, a particular subscriber unit, 

1. demodulating the signal received using a receive 
weight determined for the subscriber unit, and 
descrambling the burst payload using the CSID-based 
key for the subscriber unit; 

2. comparing the received CRC with the CRC computed 
from the demodulated descrambled bit sequence of the 
burst pay load; 

3. determining if both significantly differ, which indicates 
a transmission error or that the key is wrong, while the 
UW shows no errors, and if the condition is met, 
triggering a counter. If weight "tracking" is included, if 
the condition is not met, the communication is assumed 
not to be interferer-locked and the weight used for 
receiving from the subscriber unit (or the subscriber 
unit spatial signature) is saved ("tracked") as a "good" 
value for the subscriber unit; and 

4. if a certain number of consecutive bursts meet the 
condition stated in step 3, as determined by the counter, 
determining the port is determined to be interferer- 
locked. 

In the present PHS specification, using the unique word 
and CRC to determine interference locking would not work 
when the co-channel users are all different spatial channels 
in the same conventional channel because the CSID is the 
same for all subscriber units of the same conventional 
channel. Determining interferer-locking for the case of the 
co-channel users being the spatial channels of the same 
conventional channel can be done by maintaining spatial 
signature histories for the co-channel users. 

Once the outputs are sorted, the result is a set of output 
signals from each of the ports. 
The Apparatus 

FIG. 9 shows a block diagram of an apparatus imple- 
menting an aspect of the invention. The apparatus for 
determining the weight vector for receiving a particular 
signal transmitted by a particular subscriber unit includes 
initializing means 902 for initializing with a first initial 
vector value, first iterative means 905 for iteratively modi- 
fying the weight vector according to a first iterative method 
which minimizes a first cost function, the first iterative 
method being an iterative weight determining method with 
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good convergence properties, preferable the constant modu- 
lus method implemented as described hereinabove. The 
apparatus also includes second iterative means 907 for 
iteratively modifying the weight vector according to a 
second adaptive method which minimizes a second cost 
function, the second adaptive method being a rapidly con- 
verging iterative weight determining method, preferably the 
decision directed method described hereinabove. Initializer 
902, first iterative means 905, and second iterative means 
907 are under control of control means 911 programmed to 
activate first iterative means 905 starting from the first initial 
vector value provided by initializing means 902 until a 
switching criterion is met, the final weight vector after the 
final iteration of the first adaptive method being a second 
vector value, and to activate the second iterative means 907 
starting from the second vector value to determine the 
weight vector 909. The weight vector 909 is used by a spatial 
processor and demodulator 915 to produce, when directed 
by controller 911, the demodulated signal, the spatial pro- 
cessor using signals received via receivers 122 at the 
antenna array 103. Each of the iterative methods includes 20 
determining a copy signal. The apparatus preferably 
includes a SINR estimator 913 to estimate the post copy 
SINR in the first iterative means copy signal using the 
weight vector determined by the first iterative means 905, 
and the switching criterion preferably is the SINR estimate 
exceeding a SINR threshold 

The weight determining apparatus preferably includes at 
least one digital signal processor (DSP) devices at the base 
station, and the elements 902, 905, 907, 909, 911, 913, and 
915 preferably are implemented as programs in the one or 
more DSR 

As will be understood by those skilled in the art, the 
skilled practitioner may make many changes in the methods 
and apparatuses as described above without departing from 
the spirit and scope of the invention. For example, the 
communication station in which the method is implemented 
may use one of many protocols. In addition, several archi- 
tectures of these stations are possible. Many more variations 
are possible. The true spirit and scope of the invention 
should be limited only as set forth in the claims that follow. 

What is claimed is: 

1. A method for performance improvement of a commu- 
nication receiver receiving signals transmitted from one or 
more subscriber units, the communication receiver having 
an array of antenna elements, the method comprising smart 
antenna processing signals received by each of the antenna 
elements of the antenna array to provide a smart antenna 
processed signal, the adaptive smart antenna processing 
according to a weight vector determined from the signals 
received by each of the antenna elements, the weight vector 
determining comprising: 

initializing with a first vector value; 
until a switching criterion is met and starting from the first 
vector value, iteratively modifying the weight vector 
according to a first adaptive method which minimizes 
a first cost function, the first adaptive method being an 
iterative weight determining method with good con- 
vergence properties, the final weight vector after the 
final iteration of the first adaptive method being a 
second vector value; 
starting from the second vector value, modifying the 
weight vector according to a second adaptive method 
which minimizes a second cost function, the second 
adaptive method being a rapidly converging iterative 
weight determining method, 
each iteration of both the first and the second adaptive 
methods comprising determining a copy signal formed from 



sets of samples from each of the received signals, one set of 
samples from each received signal, the copy signal formed 
using the then current weight vector value. 

2. The method of claim 1, wherein the signals are received 
at each antenna element burst-by-burst, and wherein the sets 
of samples are from mutually contemporary bursts, and 
wherein the adaptive smart antenna processing of any set of 
mutually contemporary bursts uses the weight vector deter- 
mined from sets of samples of the same mutually contem- 
porary bursts. 

3. The method of claim 1, wherein the weight vector 
determining is blind. 

4. The method of claim 1, wherein the weight vector 
determining uses at least one digital signal processor. 

5. The method of claim 3, wherein the received signals 
comprise TDMA signals. 

6. The method of claim 5, wherein the received signals 
substantially conform to PHS signals. 

7. In a communication receiver receiving signals trans- 
mitted from one or more subscriber units, the communica- 
tion receiver having an array of antenna elements and 
adaptive smart antenna processing means, the adaptive 
smart antenna processing means including means for 
weighting the signals received by each of the elements of the 
antenna array in amplitude and phase according to a weight 
vector for each subscriber unit, the weighting forming a 
copy signal for the subscriber unit, a method of determining 
the weight vector for receiving a particular signal transmit- 
ted by a particular subscriber unit, the method comprising: 

initializing with a first initial vector value; 
until a switching criterion is met and starting from the first 
initial vector value, iteratively modifying the weight 
vector according to a first adaptive method which 
minimizes a first cost function, the first adaptive 
method being an iterative weight determining method 
with good convergence properties, the final weight 
vector after the final iteration of the first adaptive 
method being a second vector value; and 
starting from the second vector value, modifying the 
weight vector according to a second adaptive method 
which minimizes a second cost function, the second 
adaptive method being a rapidly converging iterative 
weight determining method. 

8. The method of claim 7, wherein the switching criterion 
45 is a specified first number N a of iterations, 

9. The method of claim 7, wherein the first adaptive 
method includes a copy generation step, and wherein 
switching criterion is the estimated SINR at the output of the 
copy generation. 

10. The method of claim 7, wherein the first adaptive 
method is a partial property restoral method and the second 
adaptive method is a decision directed method. 

U. The method of claim 10, wherein the first adaptive 
method is a constant modulus method. 

12. The method of claim 10, wherein each iterative 
method includes a copy generation step, and wherein the 
second cost function includes a difference term, the differ- 
ence being between a weighted signal and a decision 
directed reference signal formed from the copy signal, the 

60 decision directed reference signal forming including a track- 
ing mechanism that forms the phase of the reference signal 
at a sample point by relaxing the phase of a signal ideally 
advanced from a previous reference signal sample, towards 
the phase of the copy signal at the same sample point. 

13. The method of claim 10, wherein each iterative 
method includes a copy generation step, and wherein the 
first cost function includes a square of difference term, the 
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difference being between a weighted signal and a constant 
modulus reference signal formed from the copy signal. 

14. In a communication receiver receiving signals trans- 
mitted from a plurality of subscriber units, the communica- 
tion receiver having an array of antenna elements and 5 
adaptive smart antenna processing means, the adaptive 
smart antenna processing means including weighting means 
for weighting the signals received by each of the elements of 
the antenna array in amplitude and phase according to a 
weight vector for a particular remote subscriber unit, the 1Q 
weighting forming a copy signal for that subscriber unit, a 
method of determining the weight vectors for receiving 
signals transmitted by the plurality of subscriber units, the 
method comprising: 

for each subscriber unit, initializing with a first initial 15 

vector value, the set of first initial vector values being 

sufficiently mutually independent; and 
for each weight vector for each subscriber unit, 

until a switching criterion is met and starting from the 
first initial vector value, iteratively modifying the 2 o 
weight vector according to a first adaptive method 
which minimizes a first cost function, the first adap- 
tive method being an iterative weight determining 
method with good convergence properties, the final 
weight vector after the final iteration of the first 2 s 
adaptive method being a second vector value; and 

starting from the second vector value, modifying the 
weight vector according to a second adaptive method 
which minimizes a second cost function, the second 
adaptive method being a rapidly converging weight 30 
determining method. 

15. The method of claim 14, wherein the switching 
criterion is a specified first number N a of iterations. 

16. The method of claim 14, wherein the first adaptive 
method includes a copy generation step, and wherein the .35 
switching criterion is the estimated SINR at the output of the 
copy generation step. ^ 

17. The method of claim 14, wherein the first adaptive 
method is a partial property restoral method and the second 
adaptive method is a decision directed method. 40 

18. The method of claim 17, wherein the first adaptive 
method is a constant modulus method. 

19. The method of claim 14, wherein the set of first initial 
values are linearly independent. 

20. The method of claim 19, wherein the set of first initial 45 
values are the largest eigenvectors of R^, wherein there are 

m antenna elements and R„ is the autocorrelation matrix of 
the m-vector of signals formed at m antenna elements. 

21. In a communication receiver receiving signals trans- 
mitted from one or more subscriber units, the communica- 50 
tion receiver having an array of antenna elements and 
spatio-temporal processing means, the spatio-temporal pro- 
cessing means including means for jointly weighting and 
time-equalizing the amplitude and phase of the signals 
received by each of the elements of the antenna array 55 
according to a complex valued weight matrix for each 
subscriber unit, the convolving forming a copy signal for the 
subscriber unit, a method of determining the weight matrix 
for receiving a particular signal transmitted by a particular 
subscriber unit, the method comprising: 60 

initializing with a first initial matrix value; 

until a switching criterion is met and starting from the first 
initial matrix value, iteratively modifying the weight 
matrix according to a first adaptive method which 
minimizes a first cost function, the first adaptive being 65 
an iterative weight determining method with good 
convergence properties, the final weight matrix after 
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the final iteration of the first adaptive method being a 
second matrix value; and 

starting from the second matrix value, iteratively modi- 
fying the weight matrix according to a second adaptive 
method which minimizes a second cost function, the 
second adaptive method being a rapidly converging 
weight determining method. 

22. In a communication receiver receiving signals trans- 
mitted from one or more subscriber units, the communica- 
tion receiver comprising an array of antenna elements and 
adaptive smart antenna processing means, the adaptive 
smart antenna processing means including means for 
weighting the signals received by each of the elements of the 
antenna array in amplitude and phase according to a weight 
vector for each subscriber unit, the weighting forming a 
copy signal for the subscriber unit, an apparatus for deter- 
mining the weight vector for receiving a particular signal 
transmitted by a particular subscriber unit, the weight deter- 
mining means comprising: 

initializing means for initializing with a first initial vector 
value; 

first iterative means for iteratively modifying the weight 
vector according to a first adaptive method which 
minimizes a first cost function, the first adaptive 
method being an iterative weight determining method 
with good convergence properties, 
second iterative means for iteratively modifying the 
weight vector according to a second adaptive method 
which minimizes a second cost function, the second 
adaptive method being a rapidly converging iterative 
weight determining method, and 
a controller for 

activating the first iterative means starting from the first 
initial vector value provided by the initializing 
means until a switching criterion is met, the final 
weight vector after the final iteration of the first 
adaptive method being a second vector value, and 
activating the second iterative means starting from the 
second vector value to determine the weight vector. 

23. The apparatus of claim 22, wherein the switching 
criterion is a specified first number Nj of iterations. 

24. The apparatus of claim 22, further comprising an 
SINR estimator, wherein the first adaptive method includes 
a copy generation step, wherein the SINR estimator esti- 
mates the SINR at the output of the copy generation step in 
the first adaptive method, wherein the SINR estimator 
output is coupled to the controller, and wherein switching 
criterion is the estimated SINR at the output of the first 
adaptive method copy generation step. 

25. The apparatus of claim 22, wherein the first adaptive 
method is a partial property restoral method and the second 
adaptive method is a decision directed method. 

26. The apparatus of claim 25, wherein the first adaptive 
method is a constant modulus method. 

27. The apparatus of claim 25, wherein each iterative 
method includes a copy generation step, and wherein the 
second cost function includes a difference term, the differ- 
ence being between a weighted signal and a decision 
directed reference signal formed from the copy signal, the 
decision directed reference signal forming including a track- 
ing mechanism that forms the phase of the reference signal 
at a sample point by relaxing the phase of a signal ideally 
advanced from a previous reference signal sample, towards 
the phase of the copy signal at the same sample point. 
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28. The apparatus of claim 25, wherein each iterative 
method includes a copy generation step, and wherein the 
first cost function includes a square of difference term, the 
difference being between a weighted signal and a constant 
modulus reference signal formed from the copy signal. 5 

29. The apparatus of claim 22, wherein the signals 
received at each antenna element comprise a sequence of 
bursts, and wherein the adaptive smart antenna processing of 
any set of mutually contemporary bursts uses the weight 
vector determined from sets of samples of the same mutually 10 
contemporary bursts. 
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30. The apparatus of claim 22, wherein the weight vector 
determining is blind. 

31. The apparatus of claim 22, further comprising at least 
one digital signal processor. 

32. The apparatus of claim 29, wherein the received 
signals comprise TDMA signals. 

33. The apparatus of claim 32, wherein the received 
signals substantially conform to PHS signals. 
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[57] ABSTRACT 

A digital variable rate demodulator within a receiver oper- 
ates close to the Nyquist rate. This serves to recover correct 
timing and filter adjacent channels. The samples of an 
incoming data signal are divided into phases and combined 
into phase vectors. Intermediate points within a given phase 
vector are determined by interpolation. The data is then 
converted into a weighted sum for the purpose of decimating 
down to the baud rate. The signal-to-noise ratio is then 
optimized by estimating the likelihood of occurrence of a 
given symbol within the waveform and filtering the near 
Nyquist data rate down to a one sample per symbol data rate. 
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NEAR NYQUIST RATE VARIABLE RATE 
RECEIVER 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 5 
The present invention relates to digital signal receivers. 

More specifically, the present invention relates to a digital 
variable rate demodulator operating close to the Nyquist 
rate. 10 

2. The Background 

In a communication system, data is typically formatted 
onto a carrier signal and conveyed via a transmitter. Once 
the signal travels through some intervening medium, it is 
received and decoded by the receiver. In theory, the wave- 15 
form of the data transmitted would remain unaltered during 
the transmission process. However, in practice, the wave- 
form is distorted and corrupted by its passage through the 
electronic circuitry of the transmitter and receiver, as well 
as, through the intervening medium. Thus, the receiver 20 
serves the dual purpose of decoding and insuring that the 
received signal is free from corrupt signal data that has been 
introduced during the communication process. 

By example, in a standard satellite communication system 
a data signal is created at one location, encoded onto a radio 25 
signal, and transmitted to a satellite which may be in 
synchronous orbit above the earth. The satellite then retrans- 
mits the received signal to another location, where it is 
received and demodulated. In such a communication 
process, the data carrying signal is likely to have passed 30 
through several electronic systems, free space distances in 
excess of 40,000 miles, and twice through the atmosphere; 
subjecting the signal to numerous encounters with interfer- 
ence and distortion. 

Historically, the transmitted signal had a fixed data rate 
known to the receiver which allowed for the receiver to be 
configured for the precise characteristics of the known 
transmitted signal. Currently, the trend is toward having the 
capability to vary the data rate to accommodate the desire for ^ 
variance in satellite channel signals. For example, a single 
satellite channel may be used to carry various types of data 
signals, some of which are transmitted at low data rates and 
some of which are transmitted at high data rates. In another 
example, when the satellite signal carries compressed data, 45 
it may be desirable to vary the data rate for a given channel 
depending on the programming to be broadcasted. Thus, 
based on the differences in the speed of the action, a video 
feed of public service announcements would require a lower 
data rate than a video feed of a more data-intense activity, SQ 
such as, a sporting event. 

Current variable rate digital communication receivers, 
such as satellite receivers or cable receivers, use hybrids of 
analog and digital techniques or pure digital circuitry to 
implement the filter scheme within the demodulator. A 55 
demodulator within a digital signal receiver acts to recover 
an original signal from a modulated carrier. 

In order to provide a variable bandwidth filter structure, 
the analog/digital hybrid techniques rely on changing the 
sampling frequency and using different filters according to 60 
the required bandwidth. Thus, for a given data rate and given 
filter structure, the filter bandwidth can be adjusted by 
changing the sampling frequency. Ordinarily, this would be 
accomplished by use of a single digital match filter. This 
type of filtering scheme allows for the filter to be modified 65 
by changing the sampling frequency in order to filter vari- 
able bandwidth signals. 
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However, in using such a analog/digital hybrid scheme it 
is not viable to reduce the sampling rate below the level 
dictated by the front-end anti-aliasing filter. This level can 
also be defined as the Nyquist rate (twice the Nyquist 
frequency which is equivalent to the width of the band of 
frequencies within the waveform). In practice, when the 
sampling frequency falls below the Nyquist rate aliasing 
occurs. When aliasing occurs, frequencies greater than one- 
half the sampling rate become indistinguishable from fre- 
quencies in the fundamental bandwidth and the overall 
system is disrupted. For video communications this would 
result in unacceptable degradation to the picture and sound. 
To compensate for this anomaly at lower sampling rates, a 
lowpass filter with a lowpass bandwidth is required. This 
means that for a typical application numerous filters will be 
required to achieve the necessary output. These front end 
analog filters tend to be bulky and, thus take up critical area 
on the surface of an integrated circuit. In addition, the use of 
numerous filters in a given demodulator design readily 
becomes cost prohibitive in commercial applications where 
today's markets demand cost efficiency. 

Alternatively, systems have employed the use of pure 
digital circuitry to implement the variable bandwidth filter 
structure. The digital variable bandwidth filter structure is 
typically constructed from decimation filters. Decimation is 
a means by which the digitized video waveform is image 
scaled. By using filters to accomplish the scaling, aliasing 
concerns are limited and image artifacts are smoothed over 
allowing for further signal processing to proceed error-free. 
Filter decimation is typically accomplished by bandwidth 
limiting the image horizontally and vertically. However, 
each scaling factor requires different filter coefficients and 
adds to the complexity of the filter. Thus, in these systems 
when the need arises to tune to lower frequency rates it 
becomes imperative to use filters with narrower bandwidths. 
As the bandwidth narrows in these filters the complexity of 
the filters increases, accordingly. This results in the need to 
implement filters which are large in size and require high 
power in operation. As is the case with a analog/digital 
hybrid scheme, use of a pure digital scheme presents the 
similar problems of area consumption on a given integrated 
circuit and, ultimately, cost efficiency. The cost of imple- 
menting large complex decimation filters readily becomes 
impractical in the current commercial integrated circuit 
market. 

BRIEF DESCRIPTION OF THE INVENTION 

The present invention implements a digital variable rate 
demodulator within a receiver which operates close to the 
Nyquist rate. This serves to recover correct timing and filter 
adjacent channels. The samples of an incoming data signal 
are divided into phases and combined into phase vectors. 
Intermediate points within a given phase vector are deter- 
mined by interpolation. The data is then converted into a 
weighted sum for the purpose of decimating down to the 
baud rate. The signal-to-noise ratio is then optimized by 
estimating the likelihood of occurrence of a given symbol 
within the waveform and filtering the near Nyquist data rate 
down to a one sample per symbol data rate. 

OBJECTS AND ADVANTAGES OF THE 
INVENTION 

Accordingly, it is an object of the present invention to 
provide a filter scheme within a digital rate demodulator 
which occupies far less surface area on an integrated circuit 
in comparison to prior art filter schemes. 
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It is a further advantage of the present invention to monitor the frequency signal to approximate the DC rate, the 

provide a filter scheme within a digital rate demodulator demodulation scheme 20 and a forward error correction 

which uses proportionally less filters than prior art demodu- decoder 40 which insures that the signal, prior to channel 

lators. output 42 is free of corrupt data and matches the data output 

It is a further advantage of the present invention to 5 from me timing loop with values of symbols at the trans- 
provide a variable rate receiver which operates close to the mitter as a means of decoding. The control module 18 has a 
Nyquist rate and, therefore, because of this lower sampling first control module and a second control module, 
rate, requires far less power consumption than prior art FIG. 2 depicts a block diagram of the digital variable rate 
receivers. demodulator apparatus 20 operating close to the Nyquist 

It is a further advantage of the present invention to 10 rat e, in accordance with a presently preferred embodiment of 

provide a variable rate receiver which combines the function the present invention. A shuffler 22, a variable interpolator 

of the timing recovery with decimation filter implementa- 24 and an accumulator 26 functioning orthogonally in 

tio n . conjunction with a control module 18 implement a maxi- 

These and many other objects and advantages of the „ d ™ ted ^ sc^mc 28. The timmg recovery loop 

present invention will become apparent to those of ordinary 15 30 >s earned out by the combination of the match filter 32 

skill in the art from a consideration of the drawings and and tne "^ant interpolator 34. 

ensuing description of the invention. Focusing on the FIG. 3 block diagram illustrating the 

maximally decimated filter scheme 28 in accordance with a 

BRIEF DESCRIPTION OF THE DRAWINGS presently preferred embodiment of the present invention, the 

. 20 shuffler 22 takes the incoming signal 12 and, at each clock 

FIG. lis a system block diagram of the overall system in cycle> divides it mt0 smaller each sequence 

accordance with a presently preferred embodiment of the contaming different phases of the signal. The phases are then 

present invention. assembled into phase vectors. A shuffler 22, typically, con- 

F1G. 2 is a system block diagram depicting the digital sists of a series of levels of delays, each level having more 

variable rate demodulator apparatus in accordance with a 25 delays than the previous level. A control module 18 is 

presently preferred embodiment of the present invention. coupled to the shuffler 22 for determining the number of 

FIG. 3 is a system block diagram illustrating the maxi- samples to be taken and the distance between samples in a 

mally decimated filter scheme in accordance with a pres- given phase vector. Shufflers are well known in the prior art 

ently preferred embodiment of the present invention. and for this invention numerous designs of the shuffler can 

FIG. 4 is a waveform diagram depicting phase vectors 3 ° be employed which will satisfy the objective, 

output from the shuffler and interpolated points in between The output of the shuffler 22 is then processed through a 

in accordance with a presently preferred embodiment of the variable interpolator 24. The variable interpolator 24 in this 

present invention. scheme serves the dual functions of interpolation, as well as, 

FIG. 5 is a block diagram illustrating the timing recovery 35 decimation. The variable interpolator 24 passes all energy 

scheme in accordance with a presently preferred embodi- forward while introducing a shift in time. FIG. 4 illustrates 

ment of the present invention. a waveform diagram 50 configured from the vector data 

FIG. 6 is an illustration of the shaped waveform at the 0Ut P Ut *° m * c shm * 22 * e P 0 ?* in benveen which 

transmittal site, the waveform prior to input at the match ^ from ? c Van ^ ^tf ^polator 24 procedure. In this 

filter and waveform output after the match filter in accor- 40 ^lustniUon the solid lines 52 exempUfy the phase vectors 

dance with a presently preferred embodiment of the present ™ ip t ^ ! * ^ ** 54 

invention interpolated points in between along the waveform 50. 

„^ _'._,.„. , The determination of which points in between require 

FIG 7 is a flowchart illustrating a presently preferred m t e rpolation, known as the delay time 56, is controlled by 

method of carrying out the present invention. ^ the variab)e in t er poi a tor 24 interacting with the control 

DETAILED DESCRIPTION OF THE m ° d "If, ^ 7" D&iaK \ *?* jf^l' 

DDcmiDcn cuonniMcviT interpolator, it allows for non-constant or variable delays to 

PREFERRED EMBODIMENT by ^ the process xhe amount of ^ delay ^ 

Those of ordinary skill in the art will realize that the determined by the amount of decimation required and the 

following description of the present invention is illustrative 50 timing error. The timing error is defined as the distance 

only and not in any way limiting. Other embodiments of the between the A/D sample and the optimum sampling point 

invention will readily suggest themselves to such persons where the signal to noise ratio is minimized. Any interpo- 

skilled in the art. lator which is known in the art can be implemented in this 

The present invention implements a digital variable rate scheme, 50 lon S **> the interpolator is variable, allowing for 

demodulator scheme within a receiver which operates close 55 constructlon of filters for any delay less than, greater than or 

to the Nyquist rate. Data being input into the receiver is ec ^ al 10 ^ sam P le clock * 

modulated, non-linear shaped data which is randomly T° e phase vector data along with the interpolated data is 
encrypted in a digitized data stream. The demodulation next processed through an accumulation procedure. Here, 
process results in structured binary data with well defined the accumulator 26 adds the estimated values outputted for 
symbol boundaries. FIG. 1 depicts a system block diagram 60 a required sample and divides by the number of estimates 
of the overall system in accordance with a presently pre- taken to arrive at an averaged value for a required sample, 
ferred embodiment of the present invention. A digital signal The accumulator 26 in this scheme can be implement by use 
receiver system 10 which processes the incoming signal 12, °f a standard adder and register which are known in the prior 
includes an analog to digital (A/D) converter 14 which art- 
converts the analog signal to digital data form by sampling 65 The combination of the shuffler 22, the variable interpo- 
the signal at a specific sampling frequency, a carrier loop 16 lator 24 and the accumulator 26 working in unison through 
which interacts with a control module 18 to adjust and the directive of the control module 18 implement a maxi- 
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mally decimated filter scheme 28. The effect of which is to 
low pass data and reject ail high pass components, so that, 
a filtered waveform close to the Nyquist rate results. Such a 
minimal complexity implementation replaces the use of 
numerous decimation filters and, in doing so, frees up 5 
valuable area on the surface of the integrated circuit. 

Following the decimation filter scheme 28 the signal 
undergoes timing recovery 30 in the match filter 32 and 
constant interpolator 34 combination. FIG. 5 is a block 
diagram depicting this timing recovery scheme in accor- 10 
dance with a presently preferred embodiment of the present 
invention. Timing recovery is accomplished by adjusting the 
phase and frequency of the samples. For this invention the 
frequency and phase of the match filter 32 and the constant 
interpolator 34 are fixed by design. The function of the 15 
match filter 32 and constant interpolator 34 combination is 
to provide optimization for the signal to noise ratio at output 
and provide a measurement of the likelihood of each sym- 
bol. 

The match filter 32 which is designed according to the 20 
minimum oversampling ratio required for the highest data 
rate and which operates at the Nyquist rate for a given 
application, initiates the optimization of the signal to noise 
ratio. Since the decimation filter scheme 28 passes the whole 
spectrum, the match filter 32 serves to reshape the waveform 25 
to match the configuration of the waveform outputted by the 
shaping filter at the signal transmitting site. FIG. 6 illustrates 
the concept of how the waveform is reshaped at the match 
filter to match the waveform at the transmitting site. The first 
waveform in the sequence 62 depicts the waveform at the 30 
transmittal site, the second waveform 64 is illustrative of a 
waveform which has been transmitted, undergone A/D con- 
version and decimation filtration and the last waveform in 
the sequence 66 is the reshaped post match filter illustration. 
The constant interpolator 34 follows the match filter 32 for 35 
the purpose of introducing a time shift. The output of the 
match filter 32 is a non-integer value Nyquist rate number, 
typically 1.35 or 1.25 samples per symbol. Timing recovery 
is accomplished at the constant interpolator by synchroniz- 
ing this rate down to 1.0 sample per symbol. Thus, the 40 
constant interpolator 34 generates the optimal sampling 
points of the match filter 32 for each symbol duration. The 
output of the constant interpolator 34 results in a measure- 
ment of the likelihood of each symbol and optimization of 
the signal to noise ratio. 45 

The present invention can also be represented as a method 
for digital variable rate demodulation. FIG. 7 is a flowchart 
depicting a presently preferred method of carrying out the 
preferred embodiment of the present invention. At step 70 
the digitized incoming signal operating at or near the 50 
Nyquist rate is split into phase vectors. The phase vectors 
make up points along a waveform and in step 72 interpo- 
lation takes place in order to ascertain further points along 
the waveform in between those points represented by the 
phase vectors. In step 74 the values for a given symbol are 55 
accumulated and averaged resulting in a stream of weighted 
sums for each respective phase. The combination of steps 
72, 74 and 76 working orthogonally and simultaneously 
implement a maximally decimated filtration process. The 
signal then undergoes timing recovery in steps 76 and 78. 60 
First, in step 76 the signal undergoes a match filtration 
process to equate the waveform to the waveform transmitted 
and to predict the likelihood of a given symbol in the 
waveform and in step 78 a constant interpolation procedure 
is employed to take the rate down from Nyquist to one 65 
sample per symbol resulting in a signal which has a maxi- 
mized signal to noise ratio. 
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While embodiments and applications of this invention 
have been shown and described, it would be apparent to 
those skilled in the art that many more modifications than 
mentioned above are possible without departing from the 
inventive concepts herein. The invention, therefore, is not to 
be restricted except in the spirit of the appended claims. 

I claim: 

1. A digital variable rate demodulator comprising: 
a shuffler having an input and an output; 

a variable interpolator having an input and an output; said 

output of said shuffler coupled to said input of said 

variable interpolator, 
an accumulator having an input and an output, said output 

of said variable interpolator coupled to said input of 

said accumulator; 
a match filter having an input and an output, said output 

of said accumulator coupled to said input of said match 

filter; and 

a constant interpolator having an input and an output, said 
output of said match filter coupled to said input of said 
constant interpolator. 

2. The digital variable rate demodulator in claim 1 
wherein said shuffler further includes a second input. 

3. The digital variable rate demodulator in claim 2 further 
comprising a control module coupled to said second input of 
said shuffler. 

4. The digital variable rate demodulator in claim 1 
wherein said variable interpolator further includes a second 
input. 

5. The digital variable rate demodulator in claim 4 further 
comprising a control module coupled to said second input of 
said variable interpolator. 

6. The digital variable rate demodulator in claim 1 
wherein said accumulator further including a second input. 

7. The digital variable rate demodulator in claim 6 further 
comprising a control module coupled to said second input of 
said accumulator. 

8. A digital variable rate demodulator comprising: 

a shuffler having a first input, a second input, and an 
output; 

a variable rate interpolator having a first input, a second 
input and an output, said output of said shuffler coupled 
to said first input of said variable rate interpolator; 

an accumulator having a first input, a second input and an 
output, said output of said variable rate interpolator 
coupled to said first input of said accumulator; 

a match filter having an input and an output, said output 
of said accumulator coupled to said input of said match 
filter; 

a constant interpolator having an input and an output, said 
output of said match filter coupled to said input of said 
constant interpolator; and 

a control module with a first output, a second output, and 
a third output, said first output of said control module 
coupled to said second input of said shuffler, said 
second output of said control module coupled to said 
second input of said variable interpolator, said third 
output of said control module coupled to said second 
input of said accumulator. 

9. A digital maximally decimated filter comprising: 

a shuffler having a first input, a second input, and an 
output; 

a variable rate interpolator having a first input, a second 
input and an output, said output of said shuffler coupled 
to said first input of said variable rate interpolator; 
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an accumulator having a first input, a second input and an 
output, said output of said variable rate interpolator 
coupled to said first input of said accumulator; and 

a control module with a first output, a second output, and 
a third output, said first output of said control module 5 
coupled to said second input of said shuffler, said 
second output of said control module coupled to said 
second input of said variable interpolator, said third 
output of said control module coupled to said second 
input of said accumulator. 10 

10. A system for receiving digital signals comprising: 
an analog to digital signal converter having an input and 

an output; 

a carrier loop having a first input, a second input and an 5 
output, said output of said analog to digital converter 
coupled to said first input of said carrier loop; 

a shuffler having a first input, a second input, and an 
output, said output of said carrier loop coupled to said 
first input of said shuffler; 2 o 

a variable rate interpolator having a first input, a second 
input and an output, said output of said shuffler coupled 
to said first input of said variable rate interpolator; 

an accumulator having a first input, a second input and an 
output, said output of said variable rate interpolator 25 
coupled to said first input of said accumulator; 

a match filter having an input and an output, said output 
of said accumulator coupled to said input of said match 
filter; 

a constant interpolator having an input and an output, said 
output of said match filter coupled to said input of said 
constant interpolator; 

a first control module having an output, said output of said 
first control module coupled to said second input of 35 
said carrier loop; 

a second control module having a first output, a second 
output, and a third output, said first output of said 
second control module coupled to said second input of 
said shuffler, said second output of said second control 40 
module coupled to said second input of said variable 
interpolator, said third output of said second control 
module coupled to said second input of said accumu- 
lator; and 

a forward error correction decoder having an input and an 45 
output, said output of said constant interpolator coupled 
to said input of said forward error correction decoder. 

11. An apparatus for digital variable rate demodulation 
comprising: 

a means for dividing the samples of each clock cycle 50 
signal into a series of phases and gathering the phases 
into phase vectors; 
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a means for interpolating points in between of said phase 
vectors; 

a means for decimating down to a baud rate the data 

output from said means for interpolating; 
a means for averaging the output of said means for 

interpolating the points in between of said phase vector; 
a means for optimizing a signal to noise ratio from said 

means for averaging; 
a means for filtering a non-integer near Nyquist data rate 

down to a data rate of one sample per symbol from said 

means for optimizing; and 
a means for estimating the likelihood of occurrence of a 

given symbol within the waveform from said means for 

filtering. 

12. A method for digital variable rate demodulation com- 
prising: 

splitting a digital signal into phase vectors; 

outputting a stream of data points including interpolated 
data points in each of said phase vectors; 

decimating said stream of data points down to a baud rate; 

averaging an output of said stream of data points includ- 
ing interpolated data points; 

filtering a non-integer near Nyquist data rate down to a 
data rate of one sample per symbol; and 

estimating the likelihood of occurrence of a given symbol 
within the waveform for the purpose of optimizing 
signal to noise ratio. 

13. The method of claim 12 wherein said splitting step 
includes the sub-step of using a shuffler to split said digital 
signal into said phase vectors. 

14. The method of claim 12 wherein said outputting step 
includes the sub-steps of: 

receiving commands on which point to interpolate from a 

control module; and 
using an interpolator in conjunction with the data in each 

of said phase vectors to interpolate an estimate for the 

waveform. 

15. The method of claim 12 wherein said decimating step 
includes the sub-step of using an accumulator to decimate 
said stream of data points. 

16. The method of claim 12 wherein said averaging step 
includes the sub-step of using an accumulator to average 
said output of said stream of data points. 

17. The method of claim 12 wherein said filtering step 
includes the sub-step of using a match filter to filter said data 
rate. 

18. The method of claim 12 wherein said estimating step 
includes the sub-step of using a constant interpolator to 
estimate said likelihood of occurrence of a given symbol. 



3/24/05, EAST version: 2.0.1.4 



